Spin relaxation at the singlet-triplet crossing in a quantum dot 
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We study spin relaxation in a two-electron quantum dot in the vicinity of the singlet-triplet 
crossing. The spin relaxation occurs due to a combined efltect of the spin-orbit, Zeeman, and electron- 
phonon interactions. The singlet-triplet relaxation rates exhibit strong variations as a function of 
the singlet-triplet splitting. We show that the Coulomb interaction between the electrons has two 
competing effects on the singlet-triplet spin relaxation. One effect is to enhance the relative strength 
of spin-orbit coupling in the quantum dot, resulting in larger spin-orbit splittings and thus in a 
stronger coupling of spin to charge. The other effect is to make the charge density profiles of the 
singlet and triplet look similar to each other, thus diminishing the ability of charge environments 
to discriminate between singlet and triplet states. We thus find essentially different channels of 
singlet-triplet relaxation for the case of strong and weak Coulomb interaction. Finally, for the 
linear in momentum Dresselhaus and Rashba spin-orbit interactions, we calculate the singlet-triplet 
relaxation rates to leading order in the spin-orbit interaction, and find that they are proportional to 
the second power of the Zeeman energy, in agreement with recent experiments on triplet-to-singlet 
relaxation in quantum dots. 
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I. INTRODUCTION 

The interest to the electron spin in semiconductors 
has revived in recent years due to potential applica- 
tions in spintronicsii^ and quantum computingi^i^ On 
the way to building spin-based qubits, a number of inter- 
esting spin-related phenomena have been studied exper- 
imentally in semiconductor quantum dots. For instance, 
single-shot read-out of an individual electron spin in a 
quantum dot has been demonstrated,^'^ optically pro- 
grammable electron spin memory using arrays of self- 
assembled quantum dots has been implemented)^ en- 
tanglement of a two-electron correlated state has been 
predicted^ and accessed in transport measurements,^ 
mixing of singlet and triplet states in double quantum 
dots (due to hyperfine interaction with nuclear spins) has 
been studied , ^"'^^i^^i^'^d'^ direct access to the spin-orbit 
interaction in nanowires has been demonstratedii^ The 
use of electron spin in quantum information requires long 
spin coherence times. Identifying the mechanisms that 
govern spin decay in nanostructures is very important. 
Knowledge about the dominant mechanism can help one 
build structures with the least coupling of spin to the 
environment and therefore increase the spin lifetimes by 
orders of magnitude. 

In semiconductor quantum dotsji^ the electron spin 
interacts with the environment via a number of inter- 
actions. The spin-orbit, hyperfine, and Zeeman interac- 
tions are extremely relevant for the coherent spin dynam- 
ics. The decoherence time T2 is limited by two kinds of 
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processes: (i) spin-flip processes and (ii) dephasing pro- 
cesses. The relaxation time Ti is determined only by 
spin-flip processes. The characteristic relaxation time Ti 
is typically very long in semiconductors. Recently, Ti was 
measured in electrostatically defined lateral GaAs quan- 
tum dotSf^i^ observing Ti from 100 /xs to Is, depending 
on the strength of the external magnetic field. In self- 
assembled InGaAs quantum dotsi Ti varies from 100 /iS 
to 20 ms. Dephasing of spin occurs due to environmental 
changes that affect the phase of a coherent superposition 
of "spin up" and "spin down" states. In a Zeeman field, 
the process (i) is associated with exchanging energy with 
the environment, whereas the process (ii) is purely elas- 
tic and results from adiabatic (with respect to Zeeman 
energy) dynamics of the environment. 

There are two major mechanisms that govern the elec- 
tron spin lifetimes in quantum dots. The spin-orbit in- 
teraction is responsible for the spin relaxation (1/Ti), 
limiting the spin coherence time to T2 < 2Ti. The hy- 
perfine interaction is responsible for the spin dephasing 
and gives the main contribution to the decoherence rate 
(I/T2). The decoherence time T2 (averaged over many 
runs) has recently been accessed in experiment)^ finding 
T2* ~ 10 ns and T2 ~ 1 fis for a GaAs quantum (double-) 
dot, in agreement with theoryjiiii^ The hyperfine inter- 
action has recently been investigated experimentally in 
Refs. [miil. 

In GaAs quantum dots, the electrons are strongly con- 
fined in one direction and can be described by a two- 
dimensional Hamiltonian. The spin-orbit interaction is 
related to the absence of inversion symmetry, cither in 
the elementary crystal cell or at the heterointerface, 
and is described by the spin-orbit terms in the electron 
Hamiltonianl-'2" that are linear in the two-dimensional 
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electron momentum. These terms can be removed in 
the case of locaUzed electrons by a spin-dependent uni- 
tary transformation2iiS^i22ii2i (for the case of delocalized 
electrons see Ref. [25l ). As a result, the contribution to 
the spin-flip rate proportional to the second power of the 
spin-orbit coupling constant appears only if one takes 
into account the Zecman splitting in the electron spec- 
trum. This results in unusually low spin-flip rates. For a 
quantum dot defined in a two-dimensional electron layer, 
the linear in momentum terms in the spin-orbit interac- 
tion are dominant, provided the dot lateral size A is much 
larger than the layer thickness d (A ^ d). 

The spin-orbit interaction in quantum dots is a weak 
perturbation on top of the dot confining potential. The 
spin-orbit length Xso — the distance a bulk electron trav- 
els until its spin precesses around in the spin-orbit field — 
is typically much larger than the quantum dot lateral 
size A. The spin-orbit interaction introduces small cor- 
rections (~ A/Aso) to the spin states of the quantum 
dot and, more importantly, mediates a coupling between 
the spin and the electron orbital bath. Thus, for ex- 
ample, the electron spin couples efficiently to acous- 
tic phononsi^i^ Coupling to a functioning QPC has 
also been considered. Recent experiments^'^'i suggest 
that the spin-phonon coupling dominates the spin relax- 
ation in GaAs quantum dots in strong magnetic fields 
(B > IT). In contrast to spin relaxation, spin dephasing 
due to spin-orbit interaction is inefficient in GaAs quan- 
tum dots,— leading to T2 ~ 2Ti for a realistic system 
with A/Ago ^ 1- 

The spin-orbit interaction can also be used as a means 
of control over the electron spin in quantum dots. An os- 
cillating electric field couples to the electron spin via the 
spin-orbit interaction and produces an electric-dipole- 
induced spin resonance in the quantum dot .22. The cou- 
pling of spin to electric fields can be envisioned as being 
due to a spin-electric moment of electron<21 A dipole- 
dipole coupling between such moments is usually much 
stronger than the magnetic dipole-dipole coupling and 
occurs between distant spin- 1/2 quantum dots coupled 
with each other only through the Coulomb interactioui^? 

Two-electron quantum dots provide a unique possibil- 
ity to verify the dominant mechanisms of spin relaxation 
in quantum dots. The two-particle ground state is a sin- 
glet at zero magnetic field. The triplet level is separated 
from the singlet by a sizable energy {Ets 1 meV in 
GaAs quantum dots). In an orbital magnetic field B, the 
singlet and triplet levels can be brought to intersection 
at i? ~ 1 T in GaAs quantum dots. The singlet-triplet 
splitting, Ets = Et — -Es, was studied as a function 
of magnetic field in a number of experiments (see e.g. 
Refs. [9l l29ll30ll3lh . Furthermore, by applying an in-planc 
component of magnetic field, one can control the Zee- 
man energy Ez independently of Ets- The triplet-to- 
singlet relaxation was studied in Refs. [boIIsiIIs^IssI with 
the helpof a pulsed relaxation measurement techniquei^ 
In Ref. |30i, the spin lifetime reached up to 200 /^s and was 
limited by cotunneling processes. In Ref. [Sll. the spin life- 



time reached 2.6 ms and showed a clear E'^'^ dependence. 

In the present paper, we consider a two-electron quan- 
tum dot and calculate the triplet-to-singlet relaxation 
rates, applying the method used in Refs. [2llf2^ to the 
two-electron case. We emphasize the necessity of three 
ingredients for the spin relaxation, 

spin relaxation = spin-orbit x Zeeman x orbital bath, 

where "orbital bath" stands for any environmental bath 
that couples to the electron charge. As an example of or- 
bital bath, in this paper, we consider the electron-phonon 
interaction, which governs the spin relaxation in single- 
electron GaAs quantum dotsi^i^iii^ii^ As in the single- 
electron problem, the effect of the spin-orbit interaction 
cancels at the lowest order in the absence of the Zee- 
man field. As a result, the lowest-order singlet-triplet 
relaxation rates are proportional to E'^, as observed in 
experiment. 

In the two-electron case studied here, the Coulomb in- 
teraction between electrons plays an important role and 
the relaxation rates depend on the scale Ets^ which can 
be varied by magnetic fields. In the case of acoustic 
phonons, an enhancement of spin relaxation is expected 
at a characteristic scale Ets ^ fis/X, where s is the speed 
of sound and A is a lateral size of the quantum dot. For 
the simplicity of discussion, we do not resolve the Zee- 
man splitting of the triplet level until the end of this 
section. In the regime Ez ^ Ets ^ hs/X, we obtain a 
relaxation rate oc E^E^g in contrast to the dependence 
oc E^ found in the single-electron quantum dotSi^'i'2ii2^ 
Note that the experimental technique used in Ref. [3l| 
allows one to study spin relaxation across the singlet- 
triplet crossing and thus to probe a wide range of energy 
scales,^ < \Ets\ < 1 meV. 

The spin-orbit interaction gives rise to avoided cross- 
ings of the singlet and triplet levels in a magnetic field. 
For the linear in momentum spin-orbit interaction, we 
find that the singlet undergoes avoided crossings with 
the triplets |r±), and the energy gap which determines 
the anti-crossing is given by A ~ Ez{ri2/Xso), where ri2 
is the average distance between the electrons. ri2 is an 
increasing function of the Coulomb interaction strength 
X/a*g (ratio of dot confinement size to effective Bohr ra- 
dius), and therefore, the energy gap A is enhanced by 
Coulomb interaction. At an avoided crossing, the spin 
selection rules are maximally violated and one would 
expect a strong coupling of spins to the orbital bath. 
The anti-crossing energy A is, however, very small com- 
pared to the optimal phonon-emission energy hs/X, and 
the spin-flip rate is, therefore, extremely small at the 
avoided crossing. [The suppression occurs here mainly 
due to a small density of states of acoustic phonons at 
the phonon emission energy hujq = A.] As a result, the 
maximum of spin relaxation takes place away from the 
point of avoided crossing, at an energy Ets satisfying 
\Ets\ ~ hs/X. The spin relaxation rate has, therefore, 
a non-monotonous behavior as function of the magnetic 
field around the singlet-triplet transition. In the case con- 
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sidered here (small Zeeman energy Ez <C hs/X), a max- 
ima in the spin relaxation rate occurs on each side of the 
singlet-triplet transition, i.e. at energies Ets — ±fi-s/A. 

The Coulomb interaction has two competing effects on 
the singlet-triplet relaxation. On the one hand, it en- 
hances the mixing between the singlet and |T±)-triplet 
states. This enhancement happens, as mentioned above, 
through the renormalization of the energy gap A. On 
the other hand, the Coulomb interaction reduces the dif- 
ference between the singlet and triplet charge densities, 
|\E's(r)p and |5't(''")P- This effect is easy to understand 
on the example of a double dot. With increasing the 
strength of the Coulomb interaction the double occu- 
pancy of the electrons on one and the same quantum 
dot decreases and the charge distributions of the singlet 
and triplet states become alike. The parameter that con- 
trols the double occupancy in the double dot is the ratio 
of interdot tunnel coupling to the onsite Coulomb repul- 
sion. In the case of the single dot with two electrons, the 
role of this parameter is played by the ratio a^/X. Thus, 
in the limit a*^/X ^ 1 the charge distributions of sin- 
glet and triplet coincide. Therefore, the sensitivity of the 
singlet and triplet states to environmental fluctuations 
that couple to the electron charge is decreased by the 
Coulomb interaction. The mechanism of spin relaxation 
in the limit of strong Coulomb interaction is governed 
by virtual transitions to states outside the singlet-triplet 
sector. 

Spin relaxation in a two-electron quantum dot has been 
studied theoretically in Refs. [ssIIsgIIstIIssI . In Refs. [ssIIsgI 
the spin-orbit interaction Hamiltonian was truncated to 
the subspace of the lowest lying singlet and triplet states 
under consideration, despite the fact that the matrix ele- 
ments of the spin-orbit interaction, being proportional to 
the electron momentum and hence to the transition en- 
ergy [Pnm = [im* /h){En - Em)rnm], couplc increasingly 
strongly to states higher in energy. As a result, an impor- 
tant contribution to the spin transition amplitudes has 
been overlooked in Refs. 35,36. This contribution can- 
cels out the triplet-to-singlet relaxation amplitudes cal- 
culated in Refs. mm at the first-order of spin-orbit inter- 
action and in the absence of a Zeeman splitting. Further- 
more, Ref. m reports a magnitude of the singlet-triplet 
avoided-crossing splitting on the order of the matrix ele- 
ment of the spin-orbit interaction (A ~ hf3/X), whereas 
the magnitude of the splitting is much smaller and given 
by the product of the matrix elements of the Zeeman and 
spin-orbit interactions divided by the orbital level spac- 
ing (A ^ Ezm* (3X/h). Here, P is the coupling constant 
of the linear in momentum Dresselhaus spin-orbit inter- 
action (neglecting the Rashba interaction for simplicity) 
and we have assumed that ri2 ~ A, which corresponds 
to a weak or moderate Coulomb interaction strength. In 
Refs.H^m, the spin-relaxation rates have been computed 
numerically and the spin-orbit interaction has therefore 
been accounted for to high orders of perturbation the- 
ory in a sufficiently large subspace of states. Our re- 
sults on the relaxation rates agree with those reported 



m Refs.H^m in the regime where the second and higher 
order in spin-orbit interaction contributions to relaxation 
amplitudes can be neglected. Parametrically, this regime 
corresponds to Zeeman splittings much larger than the 
spin-orbit interaction matrix elements, Ez ^ h/S/X. 



The paper is organized as follows. In Sec. |TT1 we intro- 
duce the model Hamiltonian, which we analyze in great 
detail throughout the paper. In Sec. IIIIl we find the en- 
ergy spectrum and the wave functions of the two interact- 
ing electrons in the absence of spin-orbit interaction using 
a variational method. In Sec. IIVI we evaluate the matrix 
elements of the spin-orbit interaction with the wave func- 
tions found in Sec. IIIIl In Sec. |Vl we revisit the problem 
of energy spectrum and wave functions of two electrons, 
but this time in the presence of the spin-orbit interaction. 
In Sec. lV Al we show that without a Zeeman splitting no 
spin relaxation occurs inside any singlet-triplet subspace 
at the first order in spin-orbit interaction (i.e. second 
order for the rates) for coupling of spin to any orbital 
bath. In Sec. IV Bl we consider a finite Zeeman splitting 
and evaluate the energy spectrum and wave functions of 
the two interacting electrons in the presence of the spin- 
orbit interaction. Here, we obtain the avoided crossings 
between the singlet and triplet levels as mentioned above. 
Because of the avoided crossings, the total spin is not a 
conserved quantity and spin relaxation takes place be- 
tween the quantum dot eigenstates under the action of 
potential fluctuations (e.g. electron-phonon interaction) . 
In Sec. IV CI we consider the admixture of excited states 
to the singlet and triplet under consideration and show 
that together with the phonon potential these virtual pro- 
cesses give rise to an additional channel of spin relaxation 
inside the singlet-triplet subspace. In Sec. IVII we intro- 
duce the quantum dot relaxation rates, which we evalu- 
ate in subsequent subsections. In Sec. IVI Al we evaluate 
the singlet-triplet relaxation rates and discuss the effect 
of the Coulomb interaction on the relaxation rates in a 
variety of regimes. In Sec. IVIBI we evaluate the triplet- 
triplet relaxation rates and find that, for typical GaAs 
quantum dot with moderate strength of Coulomb inter- 
action, the rates are comparable to the spin-flip rates in 
the one-electron regime and at the same Zeeman split- 
ting. The paper contains also a number of appendices 
which complement the main text. In Appendix \^ we 
study the variational parameters introduced in Sec. IIIIl 
and used throughout the paper. In Appendix [Bl we give 
crossover formulae for important parameters governing 
the singlet-triplet crossing. In Appendix [Cl we estimate 
the fidelity of our variational method and show that it is 
exact in the limiting cases of strong and weak Coulomb 
interaction. In Appendix [D] we calculate the matrix el- 
ements of the remainder of our variational treatment. 
These matrix elements allow us to estimate the fidelity 
in Appendix [C] and can be used to improve our varia- 
tional treatment. 
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II. HAMILTONIAN 

We consider two interacting electrons in a quantum 
dot described by the Hamiltonian 



^^0 = 

Hd = 

Hso = 

Hz = 
where 



Hd + Hso + Hz, 



E 



\2m* 



K\ri - r2| 



(1) 
(2) 



^ + P^^!) + c^iP-^! - <)] , (3) 

--1,2 

• (cTi + 0-2), (4) 



1=1,2 
1 

2 



{xj , yj ) is the j-th electron radius- vector 
and pj = —ihd/dYj + (e/c)A(rj) is the kinetic momen- 
tum, with A.{xj,yj) = {^yj,Xj)Bz/2 being the vector- 
potential; (Tj — (aj, ,a^) are Pauli matrices. The spin- 
orbit coupling is given by the Hamiltonian ([3]) , where the 
term proportional to /3 originates from the bulk Dressel- 
haus spin-orbit coupling, which is due to the absence of 
inversion symmetry in the GaAs lattice; the term pro- 
portional to a represents the Rashba spin-orbit coupling, 
which can be present in quantum wells if the confining 
potential (in our case along the z-axis) is asymmetric. 
The axes x and y point along the main crystallographic 
directions in the (001) plane of GaAs. The magnetic 
field B = i3(cos(^sin6',sin93sin6',cos0) determines the 
spin quantization axis and the magnitude of the Zeeman 
splitting, Ez — gfisB, where g is the electron ^-factor in 
GaAs. The orbital effect of B is given by the component 
Bz = BcosO entering in A(r), and it allows one to con- 
trol the singlet-triplet splitting Ets = Ets{Bz)- In this 
work, we use a parabolic confining potential for the quan- 
tum dot, which is not necessarily the case in experiment. 
The use of the parabolic confinement allows us, however, 
to give numerical estimates for the spin relaxation rates 
in the presence of the Coulomb interaction between the 
electrons. We discuss the generality of our results and 
the validity of the parabolic approximation for the dot 
confinement in Sec. I VIII 

We consider the acoustic phonons as a major source of 
orbital fluctuations in the quantum dot. The potential 
of acoustic phonons reads 



C/ph(ri,r2) = 



<v 



where creates an acoustic phonon with wave vector 
q — (q||,(72), branch index j, and dispersion LOqj; pc is 
the sample density [volume is set to unity in ([5])]. The 
factor F{qz) equals unity for \qz \ <C and vanishes for 
\qz\ ^ d"^, where d is the width of the 2D layer . We take 
into account both the piezoelectric (/3qj ) and deformation 
potential (Sqj ) kinds of electron-phonon interaction.'^^ 

Apart form coupling to the acoustic phonons, the elec- 
tron spin couples also to a variety of other degrees of 



freedom in the quantum dot environment, e.g. spins of 
nuclei, excitations on the Fermi surface (in surround- 
ing metallic gates and 2DEGs), shot noise of a QPC, 
charge switches of impurities, fluctuations in the elec- 
trostatic confinement (noisy gates), _B-field fluctuations, 
magnetic impurities, etc. Different sources of spin decay 
can be relevant in different regimes and also depending 
on the type of the experimental setup. In the quantum 
dots, in which the spin relaxation has been measured so 
far,— i^iii^'^ good agreement is obtained with the theoret- 
ical predictions of Refs. [2TI123 . which are derived from the 
combined effect of the electron-phonon, spin-orbit, and 
Zeeman interactions. On the other hand, the relaxation 
rates obtained from other mechanisms, such as coupling 
to a functioning QPC,^^ nuclear spins alone^^^^^ and in 
combination with electron-phonon interaction,^'^ etc., 
show that these mechanisms can be relevant only in spe- 
cial regimes, which are usually harder to achieve in exper- 
iment. We, therefore, consider the phonons as the orbital 
bath in this paper. We note, however, that our method is 
quite general and many other fluctuating potentials can 
also be considered as long as they are Markovian with 
respect to the spin. 



III. ENERGY SPECTRUM AND WAVE 
FUNCTIONS 

We consider here the Schrodinger equation Hd\'n) = 
i?n|n), with Hd given in Eq. and find the energies 
and wave functions of two interacting electrons. We use 
the separation of variables in the Schrodinger equation 
in terms of the center of mass, R = (ri -I- r2)/2, and 
the relative motion, r = ri — r2, coordinates. The wave 
functions can then be written as |n) = \NM)\nm), with 

\NM) = RNM{R)e''^'^'^''/V2^, M = 0, ±1, ±2, (6) 

Inm) =7^„™(r)e""^'7\/2^, m = 0, ±1, ±2, (7) 

where ipi^ and ipr are the polar angles of R and r, resp. 
Note that ri r2 corresponds to ipr ^ (fr + tt. Thus, 
TO = 0,±2,±4, ... refers to singlet states (symmetric or- 
bital wave function), and m = ±1, ±3, ±5, to triplet 
states (antisymmetric orbital wave function). The radial 
component in Eq. ^ reads^ 



R 



NM 



{R) = 



2NI 



\M\ 



{N + \M\)\ A\^f\+^ ^ 



(8) 

I A'/ 1 

where is a Laguerre polynomial (N — 0,1,2,...) 

and A = ^ h/2m.*u), with uj = \/ lo^ -I- /4 and the 
cyclotron frequency = eBz/vri^c. The radial com- 
ponent in Eq. (O is given by Tlnm{r) = /„m(?')/v^, 
where the function /„m(r) obeys the Schrodinger equa- 
tion Timfnmir) = e„m/nm(f), with the Hamiltonian 



TO* 



1/4 



1 

a*gr 



*, ,2 



m UJ 



■r\ (9) 
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where a*g = h'^K/m*e^ is the Bohr radius and r = 
|ri — r2|. Another length scale in the problem is A = 
y^2h/m*uj, which gives the characteristic extension of 
fnm(r) in the absence of interaction (a|j ^ A). The 
functions fnm (r) obey the boundary conditions /„„i (0) = 
/„„i(oo) = 0, and are normalized as follows 



fn-,n{r)fn'ra{r)dr = 5„ 



(10) 



The energy spectrum of the Hamiltonian Q is given by 



E 



NMnm 



\M\ + 1) 



— [M + m) + e„ 



huj{2N - 

(11) 

where Snm are the eigenvalues of 7Y,„ in Eq. Note 
that the quantum number m enters in as a pa- 
rameter. We apply a variational method to find the 
ground state of TLm for each m. For this, we consider 
first a Hamiltonian Tim, similar to Eq. but with 



l/a^r — *■ ^/r and with uj ^ u 
the Schrodinger equation for Hr, 
set of functions 



From the solution of 
we obtain a complete 



/2r(n + 1 + ^Jm'^ + 7) r^n+y/m^+^l 
y n!r2(l + +7) 3^i+V™"+7 

xe-'^'/2^\Fi f-n, 1 + r^/}? ] , (12) 



where n = 0, 1, 2, A — ^2h/m*Cj, T{z) is the gamma 
function, and iFi(a, 6; z) is the confluent hypergeometric 
function of the first kind; 7 > and Cb are variational 
parameters. Next, we evaluate the ground state energy 
of the Hamiltonian ([H]), using fom{f) as a trial wave func- 
tion, i.e. we calculate Som = {fom\'Hm\fQm) and obtain 



£0n 



M} 1 



+ 7/2 ^ A r(l/2+ v/m2+7) 

205 r(i + v/"i2 + 7) 



Vm2+7 



-f 1 + v/w2 + 7 



-6^2 



2J)2 



(13) 



Minimizing egm in Eq. (|T3)) with respect to 7 and 
we obtain the ground state energy and wave function of 
TLm within the variational method. A detailed analysis 
of the dependence of Cj and 7 on the Coulomb interac- 
tion strength X/a*^ is given in Appendix lAl In the limit 
of vanishing Coulomb interaction {X/a*g 0), we ob- 
tained 7 = and uj = Lu, as expected. In the limit of 
strong Coulomb interaction (A/a|j ^ 1), the variational 
parameter 7 increases like 7 = (3/4)(A/2a|j)'*/'^, whereas 
Lu tends to a constant cb = •\/3w/2. 

Now we focus on the singlet-triplet crossing, between 
the singlet state lips) — |0000) and the triplet state 
\iPt) — |000, —1), which occurs with applying an orbital 
magnetic field Bz in the presence of the Coulomb in- 
teraction (the Zeeman interaction is set to zero here). 
Using Eq. ([TT|) . we obtain for the singlet-triplet splitting 



{Ets = Et — Es) 



EtS — TTf^ 



4C.2 



(14) 



where 



eBz/rn*c and S = (eqi ^ ^oo)/^- The 



Coulomb interaction enters in S = 5{X/a*g) and, for 
vanishing strength of the Coulomb interaction, we have 
S{0) ~ 1. Furthermore, in leading order, we have S — 
1 - 0FA/4a|j, for X/a*g < 1, and 6 = i(2aJj/A)2/3, for 
X/a*g ^ 1. A crossover formula and a plot for S{X/a*^) 
are given in Appendix [Bl 

Singlet-triplet splitting in two-electron quantum dots 
has been studied experimentally down to zero magnetic 
fields. ^•'^^ In relation to such measurements, we would like 
to point out the following observation. The derivative of 
Ets (considering only the ground singlet and first excited 
triplet states) with respect to lUc in the limit lOc has 
a universal value 



dE' 



TS 



dbJc 



(15) 



for all quantum dots with circular symmetry. This re- 
sult follows from the identity^ dE^/duJc = {dH/duJc)nn, 
which is valid in general for any Hamiltonian H depen- 
dent on a parameter Uc- For the two-electron case con- 
sidered here we have H — Hd 



dujc 



(16) 



1=1,2 



where we used the cylindrical gauge. Taking into account 
conservation of the 2-component of the orbital angular 
momentum for circularly symmetric quantum dots, we 
have that Ets of the ground singlet and first excited 
triplet obeys Eq. (fT5)) . 

According to Eq. (fH)) . if (5 = 1, the singlet-triplet 
crossing does not occur even for arbitrarily large Bz, 
i.e. there is no singlet-triplet crossing without Coulomb 
interaction. The singlet-triplet crossing takes place at 
Bz — Bl, due to the orbital effect of B, in the pres- 
ence of Coulomb interaction. For strong Coulomb inter- 
action, we have 5 <C 1 and the singlet-triplet crossing 
occurs at a small value of Bz, -B* <C uJom*c/e, and thus, 
one can neglect the magnetic field dependence in A and 
6. In this case, the degeneracy field is given by B* — 
25uJom*c/e, or equivalently by uj* = wo(2a|j/Ao)^/'^, 
where uj* = eB*z/m*c and Aq = yj2h/m*uia. For weak 
Coulomb interaction, the singlet-triplet degeneracy oc- 
curs at w* = 27r~-^/'^wo(2a|j/Ao)^/'^. A crossover formula 
and a plot for lo* are given in Appendix [Bl 

The accuracy of our variational method is discussed in 
Appendix [Cl We remark that our variational method is 
exact in the limiting cases of weak and strong Coulomb 
interaction. In the crossover regime {X/a*^ ^ 1), the 
overlap of the ground state wave function /om(?') with 
the exact one is close to unity (see Appendix [C| . 

Finally, the remaining part of TLm , which was not cap- 
tured by the variational method, is studied analytically 



6 



in Appendix |D] With the matrix elements calculated 
there, one can, in principle, find also the wave functions 
frtmir) of the excited states {n > 0). However, in the 
crossover regime (A/a|j ~ 1), no expansion parameter is 
present and an exact diagonalization of Hm — 'Hm + V 
is required. In this paper, we will not need knowledge of 
the excited states. 



IV. MATRIX ELEMENTS OF Hso 

The spin-orbit Hamiltonian HgQ in Eq. ([3]) connects 
different spin components and orbital wave functions. 
The matrix elements of Hso read 

(ns|iJso|n's') - * (i?„ - E„,) ^ (n|^^.|n') • (s|^,|s'), 

i=i,2 

(17) 

where and |n) = \NMnm) are, respectively, the 
eigenenergies and eigenfunctions of Hd (see Sec. lIII[) . and 
s stands for the spin quantum numbers of two electrons. 
The vector has a simple form in the coordinate frame: 
x' = {x + y)/\/2, y' = -{x ~ y)l\/2, z' = z, namely 



J -1,2, 



(18) 



where \± = h/m*{(3 ± a) are the spin-orbit lengths. 
To obtain Eq. (jl7p . we used the definition of momen- 
tum, pj = {im*/h) [Hd,rj], with Hd given in Eq. ([2]). 
Note that, due to the linear in p form of HgQ, the ma- 
trix elements (|17p vanish at the point of singlet-triplet 
degeneracy44 

The spin states of two electrons |s) are the singlet (S) 
and triplet (T) states: 



\S) 



1 



(in> 



V2 

\T+) = ITT), 

^ (ITi> 



\To) 



V2 



-UT», 

\T-) = I ii), 

HiT)), 



(19) 
(20) 
(21) 



where arrows in the first (second) place denote the spin 
state of the first (second) electron. Using the represen- 



tation in terms of R = (ri + r2)/2 and r = ri — r2, we 
write 



(22) 



3 = 1-2 



where 'S = cri + (T2 and cr = cri — (T2, and also ^ = 
(^1 +^2)72 and — —^2- Then, it is easy to see that 
the following matrix elements vanish 

{S\Hso\S) = {To\Hso\To) = {T+\Hso\T-) = 0. (23) 

Next, since the operator S is nonzero only in the triplet 
subspace and cr only between the singlet and triplets, 
all matrix elements of Hso in the triplet subspace are 
proportional to , and those connecting singlet and 
triplet are proportional to . For a magnetic field along 
Z = B/S, we thus have 



(nT±|i/so|n'r±> 
(nr±|i7so|n'To) 



= ±2i{E^~E^,){n\e\r^)-l, (24) 
= iV2iEn-En') 



x{n\e\n')-{XT^Y), 



(25) 



where X, Y and Z = I are unit vectors of a coordi- 
nate frame with Z along the spin quantization axis. One 
possible choice of X and 1^ is: X — (sin (/?', — cos (^c', 0) 
and Y = {cos ip' cos 9 , sin ip' cos 9 , — sin 6) , where ip' — 
ip — -k/A and we wrote vectors in the frame (x',y',z). 
Note that for the vector I in this frame we have I — 
(cos iy9' sin 6*, sin 93' sin 0, cos 0). The matrix elements of 
^ can be easily evaluated using the wave functions in 
Eqs. ([6]) and ([8]). To avoid extra phase factors due to the 
(7r/4)-rotation in going from (x, y, z) to {x' , y' , z), we re- 
place ipR ^ (pB/ and (pr ipr' in Eqs. ^ and ([7]), respec- 
tively. Thus, e.g., we have R± = Rx' ± iRy' = Re^^'^^t' 
and \NM) = R„™(i?)e**^v'«7\/2^. Then, the matrix 
elements {NM\^^\N'M') of the vector 



2A_ 



2A4 



-,o 



(26) 



are given in terms of the following matrix elements. 



{NM\R±\N'M') = A5m,m'±i U. 



N,N' 



\M\ + \M'\ + 1 



N 



-'iV+|M|,Af' + |A/'| 



N + N' + 1 



(27) 



r 



In particular, for the ground state N' ~ M' = 0, we have 
only (0, ±l|i?±|00) = A non-zero, which yields 



{NM\i''m - ^^Nfi5M.,±i (t*Ai\ A;1, 0) . (28) 



Next, we calculate the remaining three matrix elements 
of Hso in Eq. ([TT]). They connect the singlet and triplet 
states, and therefore, the second term in Eq. ((22)) con- 
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tributes here. We obtain 

{nTo\Hso\n'S) = i{E^- E^,) (nj^ |n') • I, (29) 

(nT±|i/5o|n'5) = T^{Er.-E^,) 

x{n\eW)■{XT^Y), (30) 

where X and Y have been defined above. Analogously 
to the previous case, we express in terms of r± = 

which allows us to evaluate (jvm\^\r\! m!) ^ provided we 
I 



know (nm|r±|n'm'). Using the wave function in Eq. ([7l), 
we obtain 

/•oo 

{nm\r±\n'm') = 6m.7n'±i rfnm{r)fn'7n'{r)dr, (32) 
Jo 

where we recalled that TZnmir) = fnmir)/^/r, with 
fnmir) given in Eq. p^ . The function fnm(i') contains 
the variational parameters 7 and cD, which are determined 
from the minimum of the ground state energy of Tim in 
Eq. ([H]), for each |m| independently. Therefore, apart 
from the explicit dependence on \m\ in Eq. I|12p . the func- 
tion /„m(r) depends on |to| also via 7 and A ~ yJl/Cj. In 
what follows, we provide 7 and A ~ -^/l/u) with an index 
m . We thus rewrite Eq. (fT2|) in the following form 



where = 1 + m? + 7™ and Am — y^^h/mF^j^. Then, evaluating the integral in Eq. (|32p . we obtaii 



45 



. /r(n + im)r(n' + im') / l + + , 2A^, 2A; 

(nm|r±|n m ) = Amm'£'m,m'±i W . ,1^./, — \ — ^ ; -n, -n ;t„,,t,n'; ~ 



, n!n'!r(tm)r(imO ^ V ^ ' ' ' ' A^ + A^, ' A?„, + A^J ' 

(34) 

where -^2(0^; /3, 7, 7'; 2;, ?/) is a hypergeometric function of two variables, see below. In Eq. (I34p . we used the notation 



Vr(im)r(tm') Ul + xi- 



(35) 



We note that the diagonal part of Amm' gives the average 
distance between the electrons in a state with n = 0, 



~ r(3/2 + v/m2+7m) 
(r)„, = Amm = Am \ , ^ (36) 



r(i + V^nM^) 



In particular, for strong Coulomb interaction, we have 
i'^)m = ^rnlm^ = A(A/2a^)i/3, in the leading asymp- 
totic order of X/ag » 1. 

The function F2{a; /3, /?'; 7, 7'; x, y) used in Eq. (f34|) is 
defined by the series^ 

where (a;)^ = r(a; -I- k)/T(x) is the Pochhammer symbol, 
and the variables x and y should obey the condition \x\ + 
\y\ < 1 to guarantee convergence. In our case, however, 
we have x = 2A^,/(A^ + A^,) and y = 2A^/(A^, A^), 
which gives — 2 > 1. Nevertheless, the series 

in Eq. (|37p contains a finite number of terms, and thus 
converges. Indeed, since the second and third arguments 



of F2{a; P, f3'; 7, 7'; x, y) in Eq. ([M)) are negative integers, 
13 = —n and (3' = — n', the Pochhammer symbols {f3)k 
and (/?'); in Eq. ([57)1 are identically zero for k > n and 
Z > n', and thus, the summation in Eq. (j37p ends at 
k = n and I = n' . 

Finally, we note that, in particular, for n = n' = Q, 
Eq. gives 



(Om|r±|Om') = Amm'^m, 



m'±l- 



(38) 



In the absence of Coulomb interaction, Eq. ([34)1 reduces 
to a simpler expression, as given in Eq. (|27p. but with 
capital letters replaced by lower-case ones. 



ELECTRON STATES WITH SPIN-ORBIT 
INTERACTION 



V. 



Using the results of Sections |IIT] and |IVl we proceed 
now to construct a perturbation theory in the vicinity 
of the singlet-triplet crossing. We assume that the spin- 
orbit interaction in the quantum dot is weak. As we show 
below, the small parameter of our perturbation theory is 
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given by the average distance between the electrons in 
the dot divided by the spin-orbit length, {r)^/Xso ^ 1; 
where {r)^ is given in Eq. p6|) and Xso is the smallest 
absolute value of X± = h/m*{Pzt a). Next, we consider 
separately the case of zero and finite Zeeman interac- 
tion. We show that, in the absence of the Zeeman in- 
teraction, the effect of the spin-orbit interaction cancels 
out at the first order of perturbation theory, similarly to 
the one-electron quantum dotS i^^i^^'^'^i^^ In the presence 
of a Zeeman splitting, we derive an effective spin-orbit 
interaction and analyze first its matrix elements between 
the states of the lowest in energy singlet-triplet subspace. 
These matrix elements lead to avoided crossings inside 
the singlet-triplet subspace and allow for the potential 
fluctuations (e.g. electron-phonon interaction) to cause 
direct transitions between the quantum dot eigenstates. 
Finally, we analyze the matrix elements, which involve 
states outside of the considered singlet-triplet subspace. 
These matrix elements allow for virtual transition to oc- 
cur, during which the spin state may be changed. We 
show that the triplet-triplet relaxation is fully governed 
by such virtual transitions. 



A. Case of zero Zeeman interaction (g — 0) 

We consider the Hamiltonian H = Hd + Hso- In 
this case, the matrix elements of the spin-orbit interac- 
tion vanish at the singlet-triplet degeneracy point, see 
Eq. p7p . and we can use perturbation theory for the 
non-degenerate case^ up to a close vicinity of the singlet- 
triplet degeneracy point. We have to make sure only that, 
within each degenerate multiplet, such as a triplet, the 
basis states are chosen correctly. For this, we write down 



the correction to the dot Hamiltonian in the second order 



1 r 
2 



(ij' 



so 



so 



where LdA — [Hd^ A\ for VA, and the z-axis is perpendic- 
ular to the plane of the 2DEG. Going to the relative co- 
ordinates in Eq. (j39p and neglecting transitions between 
levels which never intersect, we obtain 



2m*A_A 



+ 



he 



o r 



S., (40) 



where if — —id/dipn and P^, — —id/d(pr and we have 
used the cylindrical gauge A(r) — i[B x r]. Equation 
(gOl) gives a spin-orbit induced spin-splitting of the dot 
energy levels. This splitting can be viewed as an effective 
Zeeman energy which sets the quantization axis along 
the [001] crystallographic direction. The terms which we 
neglected in Eq. (^0]) are proportional to ctz = aiz — (J2z 
and thus violate the symmetry of a Zeeman interaction. 
However, since these terms are purely off-diagonal, their 
contribution goes to higher orders. 

The leading order correction to the two-electron wave 
function can be easily found using the matrix elements 
of Hso calculated in Section IIVI We choose the unit 
vectors X, Y and Z = i to point along the axes x', y' and 
z, respectively (this fixes our spin measurement frame). 
The singlet and triplet states in the first order of Hso 
then become. 



\SNMnm)' 
\T±NMnm)' 



\SNMnm) + \N M) {\n' , m - I) (|r+)a + |T_)/3) (n',m- l|r_|nm) 

n' 

+ \n', m+l) (|r_)a + \T+)l3) {n' , m + l\r+\nm)} , (41) 
\T±NMnm) ~ \SNM)^^^{\n',m ± l)a{n' ,m± l\r±\nm) + \n',m^ to =p l|rzp|nm)} 

±\Tonm)^^^"^{\N',M ±l)a{N',M ±1\R±\NM) + |7V', M T l)/3(iV', M T , (42) 



N' 



V2: 



ToNMnm)' = iToNMnm) + \nm) ^ S^UN'M + 1) (\T)a~ \T+)(3) (N' , M + 1\R+\NM) 

h ^ — ^ 



N' 



'\N'M - 1) {\T+)a- |r_)/3) {N' , M - 1\R4N M)} . 

I 



(43) 



We note that Eqs. (gl]), (|42]) and (gS]) can be rewritten 



in the following general form 
|ns>' = (l-5)|ns>. 



(44) 



S = z|^-S + -C-T = z X ^,-<T,, (45) 
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which can be viewed as a spin-dependent gauge transfor- 
mation |ns)' = exp(— 5')|ns) in leading order of S. We 
note that this transformation is identical to that used for 
a single-electron quantum dot j^^'^^i^^'^"^ with the only dif- 
ference that we sum over electrons here. Using the states 
in Eq. (|44|l . it is easy to see that the matrix elements of 
any scalar potential, such as, e.g., the electron-phonon 
interaction in Eq. ([S]), are diagonal in the spin index, 

(ns|e^t/ph(ri,r2)e-^|n's') - (C/ph)„„, <5ss', (46) 

because S in Eq. ([l5|) . being a function of coordinates 
only, commutes with scalar potentials. Thus, the spin 
degrees of freedom in the quantum dot decouple from all 
scalar potential fluctuations in the first order of Hso- 
Next, we discuss the validity of this statement for quan- 
tum dots of arbitrary shape. 

We arrived at Eqs. (I44p and (|45|) by considering the 
harmonic confining potential. For this confining potential 
and any other one possessing a center of inversion in the 
(x, ?/)-plane, the diagonal in orbit part 

5„s,ns'=^Efe)„„•('^^)ss' (47) 

i 

is identically zero for all orbital states |n), or it can be 
made so by shifting the origin of coordinates. This al- 
lowed us to choose the quantization axis for all orbital 
levels equally (along z), as required by the interaction 
(j40|) . We now consider a dot confining potential of an 
arbitrary shape, for which S'ns.ns' in Eq- (HZl) is not nec- 
essarily zero. We also allow for an arbitrary number of 
electrons in the quantum dot. As we show below, the 
transformation (|44| takes then the form 

|ns)' = (l-5s)e-^-|ns), (48) 

where S-p and Sq are some operators, respectively, di- 
agonal and off-diagonal in the Hilbert space of Hd, i.e. 
VS-p = Sv and QSq = Sq, with Q = 1 - P and 
VA — J2n ^nn\n) {n\ for yA. Here, we assume for sim- 
plicity that the Hamiltonian Hd has no orbital degen- 
eracies (otherwise one had to regroup the sets n and s 
such that En were non-degenerate). The order of the two 
transformations in Eq. ()48|) corresponds to applying first 
the non-degenerate perturbation theory (transformation 
1 — Sq) and then the degenerate one (transformation 

We start with writing down the following formally ex- 
act equality of the Schrieffer- Wolff transformation22ii& 

{Hd + Hso)e-^ ^Hd + AH, (49) 

where 5 = -S't is chosen such that Vl^H = Ai7. The 
Hamiltonian H' — Hd + AH is diagonal in the basis of 
Hd and has the same energy spectrum as the Hamiltonian 
H — Hd+Hso', the latter is diagonal in a basis related by 
exp(— 5) to the basis of Hd- Since in our case VHso = 
and thus AiJ = 0{HIq), it follows from Eq. ^ that 



[Hd, S] = Hso in leading order of Hso- In matrix form, 
the latter equation reads 

{En - En') 5„s,n's' ^^{En- i?n' ) ^ (^.Onn' ' ('^.'■)ss' ' 

(50) 

where we used Eq. (jl7|) for matrix elements of Hso- 
Clearly, for En ^ E'n, we obtain from Eq. ([50)) that 
= *Si(^j)„„' ■ ('^j)ss" "^l^ich coincides with 
what one obtains from Eq. (|45p . Further, we note that 
the fully diagonal part S'ns,ns gives rise to phase factors, 
which can be absorbed into the definition of basis states. 
These phase factors, obviously, do not affect the spin re- 
laxation, and can be chosen arbitrarily. The diagonal in 
orbit part S'ns.ns' remains undefined by Eq. (jSOp . because 
the perturbation Hso does not lift the spin degeneracy 
in the first order {VHso — 0)- Next, we split the trans- 
formation into a product, 

e-^«(l-Ss)e-^-, (51) 

where we retain only the leading order of Sq. For Sq, 
we again have [Hd, Sq] = Hso, which now gives 

Sq^{1-V)Y^i£.^-c7„ (52) 

where the factor {l — V) zeroes out the diagonal in orbit 
part. Note that Sq does not, in general, commute with 
scalar potentials and, unless S-p can be chosen such as 
in Eq. (j47|) . there can be, in principle, spin relaxation in 
the first order of Hso- 

As we have seen above, the spin splitting in the dot 
is given by Aff in Eq. (|49p and occurs in the second or 
higher orders of Hso- From Eq. one obtains at the 
leading order 

AH='^[Sq,Hso]+0{HIo)- (53) 

The operator S-p is chosen such that Aif is diagonal also 
in the spin subspacc of each orbital level. Therefore, 
Sp depends on the details of Hd, i.e. on the confining 
potential and the number of electrons in the quantum 
dot. Next we briefly discuss AiJ for the two-electron 
case considered in this paper. 

The fine structure of a two-electron quantum dot in the 
absence of orbital degeneracy (such as singlet-triplet de- 
generacy) is described by a Hamiltonian of the following 
form 

A// = A-S + ^S^,S^S„ (54) 

where A^^ and B^^ ~ B,j^ are operators diagonal in the 
orbital space of Hd- In Eq. (|54p . we did not include en- 
ergy shifts of the singlet levels, since they are negligible 
on the scale of the dot level spacing. On the same reason, 
we also neglect the triplet shifts by setting i?^^ ~ 0. 
The term A • S in Eq. (j54p can be viewed as a spin-orbit 
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induced Zeeman interaction. In leading order, A is per- 
pendicular to the 2DEG and can be found from Eq. PO]) . 
which remains valid also in this case. On average A has 
the following magnitude 

+ J = l,2 

where (...) denotes averaging over many orbital levels. 
Note that A vanishes at zero magnetic field due to time 
reversal symmetry. The remaining terms (second order 
in S) in Eq. ([51)1 describe on-site spin anisotropy. For a 
given triplet level, we can choose a coordinate frame in 
which (n|_B^i,|n) ~ (5^^, and since = 0, there re- 

main only two independent components of Bfj,^. We note 
that Bfj^iy is of higher order in Hso than A; however, 
since the spin anisotropy interaction is time-reversal in- 
variant, B^^Y?^ are the only terms left at zero magnetic 
field. 

We now discuss the relevance of the fine structure ([54|) 
to spin relaxation in the quantum dot. We have seen 
above that, for quantum dots without center of inversion 
in the {x, y)-plane, the spin-orbit interaction can, in prin- 
ciple, lead to spin coupling to potential fluctuations at 
the first order of Hso- First, however, we mention sev- 
eral cases when this coupling results in zero relaxation 
rates for two electrons and at 5 = 0. We stress that the 
statements below are formulated for the case of a linear 
in the electron momentum H so- 
il) There can be no relaxation between singlet and 
triplet states due to spin-bath coupling at the first or- 
der of Hso and in the absence of Zeeman splitting. In- 
deed, S'ns.ns' in Eq. (j47[) has no matrix elements between 
any singlet and triplet states of the quantum dot. The 
right-hand side of Eq. ((47|) can be rewritten as follows 

j = 1.2 

using the symmetry with respect to exchange of electrons. 
As a result, the singlet state is block-diagonal in coupling 
to scalar potential fluctuations [cf. Eq. (j46|l ]. i.e. a singlet 
can relax only to a singlet. 

(ii) There can be no relaxation between spin states be- 
longing to one and the same orbital level in the first or- 
der of Hso without Zeeman interaction. This is obvious 
since P[S'g, J7ph] = 0, see Eq. ^ with S ^ Sq. 

After (i) and (ii) we are now left only with transitions 
between triplet states of different orbital levels. Let the 
triplet states be fine-split with some characteristic energy 
Afs. This fine splitting can be calculated from the Hamil- 
tonian (j54p . using coupling constants obtained from the 
perturbation theory expansion of Eq. (|49)) . However, for 
the sake of making a simple argument here, we assume 
Afs to be an independent parameter. Clearly, if Afs = 0, 
we are free to construct any linear combination of states 
within a triplet; so we can recover Eqs. ((44l) . (|45p and 
(I46p . which yield no spin relaxation. To be more rigor- 
ous, we formulate an additional condition when there is 
no spin decay in the first order of Hso- 



(iii) If the orbital relaxation time r is much shorter 
than h/ Afs, then the spin relaxation between triplet 
states belonging to different orbital levels is suppressed. 
As mentioned above, note that, for Afg <^ fi/T, one can 
choose to work with the states in Eqs. (|44)) and ([45]) 
instead of the true eigenstates pS)) . because the inter- 
nal evolution of the states occurring due to the 
fine structure splitting Afs, is slower than the lifetime 
T. Thus, the spins can be considered as decoupled from 
the potential fluctuations during the time ~ ft/Afs. As a 
result, the spin relaxation occurs due to coupling at the 
second or higher orders of Hso- 

Finally, we conclude that spin relaxation due to spin- 
bath couplings at the first order of Hso and in the ab- 
sence of the Zeeman interaction is possible only between 
spin triplet states belonging to different orbital levels and 
only if Afs > H/t. The interaction which causes spin re- 
laxation is proportional to {^^^ — ^mm) ' where n and 
m denote the orbital states involved in the relaxation 
process. We do not consider this mechanism in further 
detail here, since it refers to an asymmetric quantum dot, 
which goes beyond our model. 

B. Case of finite Zeeman interaction (g ^ 0): 
Degenerate perturbation theory 

We consider now the Hamiltonian H = Hd + Hz + Hso 
in the vicinity of a singlet-triplet crossing. We assume 
the Zeeman spitting Ez — gi^-sB to be large compared 
to both the fine structure splitting Afs and the level 
broadening H/t. This allows us to set the spin quan- 
tization axis along the applied magnetic field B. In Sec- 
tion |V3 we found the wave functions of two electrons 
in the quantum dot for the case without Zeeman interac- 
tion. It is convenient now to use these functions as basis 
states for studying the effect of Zeeman interaction. We 
perform the following unitary transformationf2iiS^i2^i2i 
H ^ e^'2{Hd + Hz + Hso)e'^'^, with Sq given in lead- 
ing order in Eq. (|52p . After this transformation, a basis 
state |ns), associated with the Hamiltonian H, will corre- 
spond to the basis state |ns)' = (1 — S'g )|ns), associated 
with the Hamiltonian H = Hd + Hz + Hso- In first order 
of Hso-: it is straightforward to obtain 

H = Hd + Hz+Hl", (57) 

^1° = EzJ2[lxQi^]-a,, (58) 

where I = B/_B is the spin quantization direction, and 
Ez = gfJ-sB is the Zeeman energy. The projector Q 
zeroes out the diagonal part of ^j. Below, we consider 
the harmonic confining potential, for which we achieve 
Q^j = by choosing the origin of coordinates in the 
dot center. In Eq. (j58p . we can rewrite the sum over two 
electrons in terms of the relative coordinates, 

J2[l>^^,]-'r,^[lxe]-^ + l[l>^e]-^- (59) 
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The first term on the right-hand side of Eq. ([59)) admixes 
only triplet states with different orbital wave functions 
and is, therefore, responsible for triplet-to-triplet relax- 
ation. This term can be taken into account by means of 
perturbation theory for the non-degenerate case. In con- 
trast, the second term on the right-hand side of Eq. ((59|) 
is strong at the singlet-triplet crossing, since it connects 
singlet and triplet states which are degenerate. For this 
term, we use perturbation theory for the degenerate case. 
The matrix elements of H§'-' in Eq. (|58p. therefore, read 



(nT±|ij|0|n'ro> 
{nT±\Hl^\n'S) 
{nTo\Hl^\n'S) 



TtV2Ez{n\^''\n') -{Xt^Y), (60) 



V2 



Ez{n\e\n')-{XT^Y), (61) 



(nr_|ff|0|n'T+ 



0, 



(62) 



where X and Y are unit vectors in the plane perpendic- 
ular to the applied magnetic field B (see Sec. IIV[) . Note 
that, according to Eq. ([S^ . there is no admixture be- 
tween any singlet \S) and triplet \Tq) states, as well as 
between any triplet states: |r+) and |T_). This will have 
an effect on the spin relaxation rates in Sec. IVI Al where 
the corresponding rates for the transitions 5 ^ Tq and 
r_ <-> T_|_ are found to be zero in leading order in spin- 
orbit interaction. Next, we use Eqs. (|60p and (pTjl to find 
the two-electron wave functions for the lowest in energy 
singlet-triplet crossing in the quantum dot. 

We consider first a close vicinity of the singlet- 
triplet crossing, where the contribution of strong ma- 
trix elements is dominant. Substituting Eq. (|3ip into 
Eqs. (ISni)-(IS21) and using Eq. §Ei), we obtain for |*s) = 
lOOOOS*) and |*To,±) = |000, -1, To,±) the following ex- 
pression 



±- 



1 

a7 



+i 



(vI/5.T±|i?|°|*To) = (vl/Tji?n«'T_) =0, 



1 COS 9 
-T^)sm^- 

SO 



cos if 

(63) 
(64) 



where we used the convention for X and Y given below 
Eq. (ESI). Next, introducing A± = | (^fT^ |ij|'^|^'s) | , we 
obtain 



Ah 



|-Ez|Ai,o 
2V2 



A 



2L 



(65) 



where I = (cos sin 6*, sin sin 0, cos 0) gives the de- 
pendence of A± on the direction of the magnetic field. 
Note that, since we are considering a particular sin- 
glet and triplet, namely |*s) = |00005) and |*To.±) = 
|000, — 1, To,±), the directions of I, for which these lev- 
els intersect each other in the absence of Zeeman and 
spin-orbit interactions, correspond to I lying in the north- 
ern hemisphere {0 < 7r/2). For I in the southern hemi- 
sphere (9 > it/2), a different triplet, |000, 1, ro_±), lies 
lower in energy than the considered one. Nevertheless, 
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FIG. 1: Sketch of the singlet-triplet crossing in a quantum dot 
witli spin-orbit interaction. The triplet level splits in three 
due to the Zeeman interaction. The singlet |^'s) undergoes 
avoided crossing with the triplets |\E't-^), with the splittings 
A± given in Eq. (|65|l . At the same time, the degeneracy of 
I'^s) and I'^Tq) at the crossing point is not lifted. 



we shall follow the triplet |*To.±) = |000, -1, To.i) to all 
values of B and directions of I, since it makes the fol- 
lowing consideration simpler. Note that the results for 
the other triplet, |000, 1, To,±), are obtained by replacing 
Bz —Bz, which, in particular, corresponds to replac- 
ing A± ^ Azp. 

The quantities A± in Eq. (p5|) give the strength of 
the singlct-triplct mixing. The Coulomb interaction en- 
ters in Eq. ([63]) through the length scale Ai.q. Note that 
Ai,o ~ Ao,Oj''^i,i for all values of A/a|j. One can inter- 
pret Ai_o as an average distance between the electrons, 
Ai.o — (r). In the limit of strong Coulomb repulsion, 
all three quantities Ao,Oj Xi^, and A14, converge to one 
another and become equal to (r) = \{\/2a'^y^^. The 
average distance between the electrons is increased by 
the presence of the Coulomb interaction (see Eq. ([S^ V 
Therefore, we expect that the Coulomb interaction gen- 
erally enhances the effects related to the spin-orbit inter- 
action. In particular, this enhancement is expected to be 
strong in quantum dots with large X/a*g. 

The singlet-triplet crossing is shown schematically in 
Fig. [TJ The triplet splits into three levels, i?T± = Et T 
\Ez\ and Etq — Et, due to the Zeeman interaction. We 
assume a negative ^-factor, i.e. Ez = — \Ez\, which is the 
case for GaAs. The avoided crossing of the singlet state 
l^^s) with the triplets \'i^T±) occurs due to the interaction 
H§'-^ , with the matrix elements given in Eq. (p5)) . The 
splitting energies are the doubled A± in Eq. ([55]) . Note 
that the degeneracy of the singlet state with the triplet 
I 'J To) at the crossing point is not hfted. This degeneracy, 
however, can be lifted in higher orders of Hso- Next, we 
consider the interaction between the states: \'^s), I^t+), 
and |5't_)- And since the state I^'tq) couples to neither 
of these states, see Eq. ([5^ and the discussion thereafter, 
we disregard it here. Writing the wave function in the 
form 



0^) 



(66) 
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we obtain the following set of equations 

Es — E WsT^ 
Wt+s 
Wt_s 




0, (67) 



where Wst± stands for s\Hz^\^t±) ■ The character- 
istic equation reads 



Eq — E - 



0, 



(68) 



where i5T± = Et T \Ez\, and Es and Et are the ener- 
gies of singlet \^ s) and triplet state I^'tq), respectively. 
Since A± ^ Ez, we can solve Eq. ([55]) in the secular 
approximation. Setting in turn A_ and A+ to zero in 
Eq. ([68]), we find the following expressions, for the lower 
(-I-) and upper (— ) solid curves in Fig. [1] 



E^ 



E' 



Est 



E 



Es) 



4A| 



(69) 



Next, multiplying Eq. §SI by {Et+ - E){Et_ - E) and 
dividing the obtained polynomial by E^ — {E+ +E^)E + 
E+E-, we obtain in leading order the solution for the 
S-shaped curve in Fig. [l] 



En 



Ex 
1 



{Et+ - Esf + 4A2. 



-^I{Et_-EsY 



-4A2 



(70) 



We note that, in the case of positive g-factor, i.e. Ez — 
\Ez\, one should modify Eqs. ([55)1 and ([70)1 as follows: 
Et± — > Et^ and A± Aip. Finally, the wave function 
of a level with the eigenenergy E = E± , Eq is given by 
Eq. ([5^ . with the coefficients 



E — Ex_ 



-a, 



\Hlo 



1 



+ 

Al 



E 



Et- 

1/2 



-a. 



{E - Et+ 



{E-Et_Y 



(71) 



Equation ((66)) together with Eqs. (|69|) — (|7T|) allows one 
to calculate transition rates between quantum dot levels 
under the action of a perturbation. As a perturbation, 
in the next Section, we consider the electron-phonon in- 
teraction. 

Now, we return to Eq. (j65p and consider the par- 
ticular case when the magnetic field B is perpendic- 
ular to the 2DEG plane [6 — 0). Then, according 
to Eq. (inSl), we have A+ = \Ez\Xi,Qm*\a\lh^/2 and 
A_ = \Ez\Xifim*\l3\/h^/2. These expressions, in prin- 
ciple, allow one to access the Rashba (Dresselhaus) cou- 
pling constant a {(3) separately in a measurement of A+ 
(A_). Accessing a and (3 separately was, so far, possi- 
ble only by analyzing the weak (anti)localization data in 
diffusive 2DEGs.'*^ The splittings 2A± can, in principle. 
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FIG. 2: (Color online) Angular dependence of the splitting 
A+ as given by Eq. (|65|) . evaluated for A+/A_ = 2. The 
direction in space corresponds to the spin-quantization axis, 
given by Z = B/B. At the "magic angles", the value of A+ 
equals zero. The angular dependence of the splitting A_ in 
Eq. (|65p is obtained from this figure by refiecting it in the 
(a;', j/')-plane. Note that the circular symmetry of the quan- 
tum dot was essential for obtaining this angular dependence. 



be deduced from transport spectroscopy of the quantum 
dot. GaAs quantum dots have, however, a fairly small 
splitting. To give an estimate, we take typical GaAs 
values: Ai.o = 100 nm for the inter-electron distance, 
h/m*f] = 8 /xm for the spin-orbit length^ and a mag- 
netic field _B = lOT (with g — —0.44 ); as a result, we 
obtain 2A_ « 20 ^eV. A different possibility to access 
A± is to measure the spin relaxation times and to de- 
duce the strength of the spin-orbit interaction therefrom. 
We also note that an analogous spin splitting (though 
of a much larger value) has been observed by means of 
transport spectroscopy in a recent experiment on an elon- 
gated quantum dot defined electrostatically inside a InAs 
nanowire4^ 

Equation ([^ has a strong dependence on the direc- 
tion of Z = B/B. For some angles — "magic angles" 
— the effect of the spin-orbit interaction vanishes at the 
leading order. From Eq. (p5|) . we find that A+ = for 
cos 9 — A_/A+ and ip' = 0, tt, provided |A+| > |A_|. 
Alternatively, if |A-| > |A+|, we find that A+ = for 
cos 61 = A+ZA- and (p' — 7r/2,37r/2. Similarly for A_, we 
have the same expressions, but with cos 9 ^ — cos 9. In 
Fig. [2l we plot the angular dependence of A_|_ given in 
Eq. (|65p. calculated for the ratio A+/A_ = 2. A similar 
figure can be obtained also for A_ in Eq. (p5)) . by reflect- 
ing Fig.[2]in the (x, ?/)-plane. Clearly, the "magic angles" 
can be used to determine the relative strength between 
the Rashba {a) and Dresselhaus (/3) coupling constants 
(recafl that 1/A± = m*(/3 ± a)/h). 

The angular dependence of the splittings Aj- in 
Eq. (|65p is specific to the circular symmetry of the har- 
monic confining potential used above. In the opposite 
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FIG. 3: (Color online) Opposite to Fig. [21 we plot here the 
angular dependence of the splittings A± (with A+ = A_) as 
given by Eq. (|72|l . which describes a quantum dot without 
circular symmetry. The degree of asymmetry needs to be 
large enough, such that in a magnetic field causing the singlet- 
triplet crossing the imaginary part of f = {tpT\F\tps) to be 
negligible. For the plot we assumed fyi\+/f^iX- = 2. 



extreme case, when the two-electron wave functions can 
be chosen real (or their imaginary parts due to Bz are 
negligible) we obtain A_|- = A_ and 

where f = {iI:t\^\'4^ s) is the dipolar matrix element (not 
to be confused with (r) ) , which we assume to be real. We 
note that this corresponds to a quantum dot elongated in 
one direction, or to a double quantum dot. The singlet- 
triplet crossing can occur in this case at a value of 
much smaller than the characteristic confining magnetic 
field. Then, the imaginary part of f is negligible. Note 
that Eq. ([72]) can be rewritten in the form 



A± = ^W^ + 5|ly'l - cos2(^' - <^o) sin^ 0, 



with iff) satisfying tan Lp^ — fx'\- /fyiX-^-. Then, it is 
easy to see that Eq. ([72]) describes a doughnut (or more 
precisely a horn torus) with the rotation axis at Lp' — ipo 
and 9 = 7r/2, see Fig. [31 In this case, the "magic angle" 
coincides with the rotation axis of the torus. Knowing 
the rotation axis, one can find the relative strength of a 
and /3, provided the asymmetry of the quantum dot in the 
directions x' and y' is known. In Fig. [31 we plotted the 
angular dependence of A± assuming fyi\j^jfx'\- = 2. 



C. Case of finite Zeeman interaction (3 7^ 0): 
Non-degenerate perturbation theory 

In Sec. IV B[ we have taken into account the matrix 
elements (|63l) by means of the degenerate perturbation 
theory and we obtained two avoided crossings between 
the singlet and the triplet |4't±)- We expect that ac- 
counting for these (strong) matrix elements should be 
sufficient in the majority of cases and it should give the 
dominant channel of spin relaxation in the vicinity of the 
singlet-triplet crossing. To complete our analysis of the 
first order in Hso contributions to the spin relaxation, 
we must also consider the rest of the matrix elements 
in Eqs. ([SO)) and ([CT|) . We call these matrix elements 
weak. The weak matrix elements, despite the fact that 
they do not induce strong mixing between singlet and 
triplet states, can, in principle, govern the spin relax- 
ation in special cases, when the strong matrix elements 
become, on some reason, inefficient. Examples of such 
cases include: positions far away from the singlet-triplet 
crossing, where Ets is comparable to the energy distance 
to other excited states; magnetic field B along a "magic 
angle" ; and the limit of strong Coulomb interaction. 

We use perturbation theory for the non-degenerate 
case and account for the admixture of distant in energy 
states to the singlet-triplet subspace under consideration. 
For these we first separate the weak matrix elements from 
the Hamiltonian into a perturbation (cf. Eq. ([57]) ). 

H = Ho + V^, (73) 
Ho - Hd + Hz+VsTHl^, (74) 
Vw = {1-Vst)H^z°, (75) 

where Vst projects onto the singlet-triplet subspace, 
VstA = iZ^u^t^AtAivl with (*s,*T±,*To) and 

VA. Note that all the strong matrix elements are in- 
cluded into Hq, whereas all the weak ones into V^j- Next, 
we transform H to H' = (l + 5w)-ff(l-5'w), with 5w < 1 
obeying [Hq, S^,] = Kv- After this transformation, we ob- 
tain H' = Hq + 0{HgQ). At the same time, by doing 
this transformation, we correct the basis states to take 
into account the perturbation V„ to leading order. 

Considering now the phonon potential C/ph(ri, r2, 
we apply the same unitary transformations, C/pj^ = (1 -I- 
'S'w)C^ph(l — '5'w), and obtain 

C/ph^C/ph + [^w,C/ph]. (76) 

The transformation matrix can, formally, be written 
as follows 

S^ = izJ^Hl^ + 0{Hlo), (77) 

where L^A = [i/d,^] and LzA = [i/z, A], for VA. For 
simplicity, we assume small Zeeman splittings, Ez 
hwQ^ which is usually the case for GaAs quantum dots. 
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Neglecting then the Liouvillean Lz in the denominator 
of Eq. ([77)1 and going to relative coordinates, we obtain 

= Ez X . S + X C\ ■ rX . 

(78) 

The first term in Eq. ((78|) is easy to evaluate for the har- 
monic confinement and because the Coulomb interaction 
does not couple to the coordinate R. It gives 

|i[Zx|«].S = xD«].S + /(R), (79) 

:= {\Z^dldRy,,X-^^dldR^,,Q) , (80) 

where /(R) is a function of coordinates only, which drops 
out upon substitution of into Eq. ((76|) . The second 
term in Eq. (|78p is more difficult to evaluate and it is the 
subject of a separate publication. 

We note here that the first term in Eq. ([75]) is responsi- 
ble for the triplet-triplet relaxation. In combination with 
the phonon potential (or any other orbital bath), it re- 
sults in finite relaxation rates for the transitions r+ ^ Tq 
and T_ <-> Tq. The transition T+ <-^- T_ is, however, for- 
bidden at the first order of HsOi because it requires a 
change of spin AS^ — 2. We evaluate the triplet-triplet 
relaxation rates in Sec. IVI Bl 

The second term in Eq. (|78p is similar by nature to 
the first term, however, it gives rise to singlet-triplet re- 
laxation. The projector 1 — Vst^ multiplying the sec- 
ond term, ensures that only distant in energy states are 
admixed to the singlet-triplet subspace under consider- 
ation. In the vicinity of the singlet-triplet (avoided) 
crossing, there is a strong admixture between the sin- 
glet and triplet states inside the singlet-triplet subspace, 
see Sec. IV Bl Therefore, the second term can usually be 
neglected, except for special cases mentioned above. The 
second term contributes, due to its spin symmetry, to 
the transitions Tj^ <-!■ S and T_ <-> S, as also the strong 
matrix elements do. In this paper, we neglect the contri- 
bution of the second term to the spin relaxation. 

VI. PHONON-INDUCED SPIN RELAXATION 

The spin dynamics, decoherence, and relaxation are 
described by the quantum dot density matrix p(t). 
Within the Bloch-Wangsness-Redfield theory;^ the den- 
sity matrix obeys a set of linear differential equations, 

-TTPaa'it) = -iuJaa' Paa' (t) + Rgg'fjp' Pfifj' jt) , (81) 

f)f3' 

where Raa'fSfi' = {Ra'ai3'f3)* are time- independent con- 
stants, which form the Redfield tensor, and oJaa' = {Ea — 
Ear)/h are the transition frequencies, with Ea being the 
energy of the level a. Equations ([8T|) are usually derived 
microscopically using the Born-Markov approximation.^^ 
In this approximation, the system-bath coupling is as- 
sumed to be sufficiently weak, so that the next-to-leading 



order Born terms contribute at times when the density 
matrix has already reached its equilibrium value within 
the desired accuracy. Additionally, the bath correlation 
time Tc is assumed to be much shorter than the character- 
istic times occurring in Eq. (|8T|) . As a matter-of-fact, the 
evolution of the density matrix on times shorter or com- 
parable to Tc is not described by Eq. (|8ip and is rarely 
of practical interest. For the phonon bath considered in 
this paper, the correlation time is extremely short. It 
is the time it takes a phonon to traverse a distance equal 
to the quantum dot size, i.e. Tc ~ {t)/s, where (r) is the 
average distance between the electrons and s is the speed 
of sound in the sample. For GaAs quantum dots, we es- 
timate Tc ~ 10 ps. At the same time, the characteristic 
spin relaxation times in GaAs quantum dots range from 
100 /iS to 1 s, as measured experimentally.^'^'^''^*''^^''^'^ 

In the Born-Markov approximation, the Redfield ten- 
sor reads 

Raa'pP' = '^%a'al3 + 13' a' aP 

7 7 

where (r+^,^^,)* = F^,^^,^ are given as follows 

1"°° (jf 

n^pp'^ ^e— ^^'*(a|C/i„t(OI«')(/3|C/int|/3'), (83) 



r-„./3,'= «"'*(a|C/i.t|a')(/3|t^int(OI/3')- 

(84) 

Here, the bar denotes averaging over the bath degrees of 
freedom, |a) is a state of the quantum dot, and U{at is 
the interaction between the quantum dot and the bath. 
The time dependence of Uint{t) is taken with respect to 
the bath only, C/int(0 — exp{iHBt / h)Uint eyip{—iHBt / h) , 
where Hb is the Hamiltonian of the b ath. In the expres- 
sions above it is assumed that Uint{t) — 0. 

In our case, the system-bath interaction is obtained 
from the transformation in Eq. (I76|) 

Uintit) = Uph{t) + [5w, C/ph(0], (85) 

with 5w given in Eq. (|78p . The time dependence of the 
phonon potential Uphit) is governed by the Hamiltonian 
of the free phonons, 

Hb = i/ph = ^<ii (^L^^ + 1/2) , (86) 

where the index j runs over the three branches of acoustic 
phonons. The optical phonons can be neglected, since 
usually |£^rs|, \Ez\ < ^opt, where fiwopt is the energy of 
optical phonons. The acoustic phonon potential can be 



rewritten in terms of relative coordinates as follows 

[/ph(R,r,i) = ^A/q,F(g,)e^iii-^cos(q|| •r/2) 
qj 



- \ (e/3q, -zgSq,), (87) 

y Pc^qj 

where bqj(t) = &qj exp(— iwqjt) gives the time depen- 
dence mentioned above for [/ph(i)- The averaging over 
the bath (denoted by bar in Eqs. ((83|) and (I84p ) is eval- 
uated as a thermodynamic mean using the free phonon 
statistical operator 

Pph = Z^^e-"^-/^, (88) 

where T is the temperature in units of energy (ks = 1) 
and is the phonon partition function. 

In addition to the Born-Markov approximation, a fur- 
ther important simplification arises in the case of the 
spin-phonon coupling considered here. The elements of 
the Redfield tensor Raa'pp' are much smaller than the 
corresponding transition frequencies LOaa' for a ^ a' . 
One can then expand in terms of Raa'pp' l^aa' ^ 1 in the 
characteristic equation of the master equation ([5T|) . At 
the zeroth order, the off-diagonal part of p(t) decouples 
from the master equation (jSip and can be considered sep- 
arately. The diagonal part, then, obeys the Pauli master 
equation 

= ^ [ra/3/0/3/3 — TfjaPaa] , (89) 

where Tap is obtained from Eqs. ([83]) and ((84)) as follows 

In our notation. Tap is the probability per unit time to 
transit from state /3 to state a. For the rates Tap the 
time integral can be extended to — oo in the lower bound, 
arriving at the Fermi "golden rule" rate 

^o.p = J ^e-*"°''*(/3|C/int(i)l«>(«|C^int|/?). (91) 

Provided the bath is at thermal equilibrium (which is 
our assumption), the rates Fq,^ obey the detailed balance 
equation 

Fa;3^F^^exp(^ ^^~^° ^ . (92) 

In what follows, we calculate the relaxation rates for 
transitions between the quantum dot eigenstates, using 
the system-bath interaction ([55]) and the definition ([91]) . 



15 




FIG. 4: (Color online) Spin relaxation processes accounted for 
(a) and neglected (b) in our treatment. The arrows show tran- 
sitions between states of the lowest in energy singlet-triplet 
subspace. The excited states (upper row) are involved only 
in virtual transitions as intermediate states, (a) We take into 
account the strong matrix elements, which mix the singlet- 
triplet subspace inside itself, allowing phonon-induced transi- 
tions to occur between the singlet and two triplet states. We 
denote these transitions by the lower-most arrows, where the 
product H§^Uph next to the arrows denotes that the transi- 
tion occurs under the simultaneous action of H§'^ and t/ph. 
Additionally, we take into account the weak matrix elements 
for the triplet-triplet transition. In contrast to the singlet- 
triplet transition, a triplet-triplet transition requires the par- 
ticipation of an excited state. During the virtual process, 
H§'^ and Uph act separately on the two halves of the process 
(e.g. going up with H§'^ and coming down with Uph)- We de- 
note this difference by placing and Uph on opposite sides 
of the arrow. Note that H§^ connects states with different 
spin, whereas Uph with the same, (b) Transitions between the 
singlet and triplet states due to the weak matrix elements are 
neglected. These processes are relevant only in special cases 
mentioned in the text (see Sec. IV Cp . 



Before proceeding, we summarize our main results from 
Sec. |V] regarding the admixture processes contributing 
to spin relaxation. The admixture processes are illus- 
trated in Fig.m where the arrows denote finite relaxation 
rates occurring either due to the direct mixing inside 
the singlet-triplet subspace (Sec. IV B[) or due to virtual 
excitations to states outside the singlet-triplet subspace 
fSec. lVC]) . In Fig.|4Ka), we show the processes which we 
take into account in our treatment, whereas in Fig. SJb) 
we show the processes which we neglect. The neglected 
processes were discussed in Sec. IV CI Finally, we remark 
that the following relaxation rates are identically zero at 
the leading order of spin-orbit interaction 

TsTo = rT+T_ = 0, (93) 

as well as the reverse rates which follow from the detailed 
balance 



A. Singlet-triplet relaxation rates 

The Hilbert space of our system is spanned by = 4 
states: jvPs), |*t+), |*t_) and \^To)- However, all of 
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these states, except |^'to)i are mixed with one another 
by the interaction in Eq. In Sec. IVBl we 

found the eigenstates and the corresponding energies of 
the quantum dot in the presence of H§'^. The eigenstates 
are given in Eq. 166]) . with the coefhcients in Eq. ((7T|) . 
and the renormahzed energy levels are given in Eqs. ((69|) 
and (|70p . Here, we calculate the relaxation rates between 
these eigenstates, taking into account only the mixing 
due to the strong matrix elements of . 

We proceed by letting U\at = C^ph and focusing on the 
states 



l*o), |*+>, 1*- 



(94) 



which are linear combinations of the singlet \^ s) and 
two triplet |^'t±) states, as shown in Eq. Then, it 

is straightforward to obtain the matrix elements of ?7ph, 



(^-clC/phl^-/?) 



a-aap [{'>i^s\Uph\ips) 

+ '5a/3(V'T|C^ph|'0T): 



(V'TlC/phlV'T)] 

(95) 



where a = 0, ± labels the considered eigenstates, and 
l'0s(T)) is the singlet (triplet) orbital wave function. The 
coefficients a„ read 



-1-1/2 



{E,,~Et^Y {Eo,~ET_f 



(96) 



where E^ — Eo,E± are the exact solutions of Eq. (p5)) . 
The last term in Eq. (|95p is irrelevant for our further dis- 
cussion, since we consider only Tap with a ^ (3. Thus, 
from the first line in Eq. ((95l) . we see that this mecha- 
nism of spin relaxation relies on the difference between 
the singlet and triplet charge profiles, |i/'s(ri, r2)P — 
|V'T(ri, r2)p. We show below that the Coulomb inter- 
action has a strong effect on the magnitude of this differ- 
ence, making the singlet and triplet charge distributions 
look the same in the limit of strong Coulomb interac- 
tion (A/a|j S> 1). Thus, we expect this mechanism to 
be most efficient in quantum dots with weak and mod- 
erate Coulomb interaction {X/a*g < 1), and less efficient 
in quantum dots with strong Coulomb interaction and in 
double quantum dots. 

Using the wave functions found in Scc. lIIIl we evaluate 
next the following form-factor 

^ST(q||) = (^5|e'''«^cos(q||r/2) I7/-5) 

-(^T|e^'i|i^cos(q||r/2)|^T). (97) 

For the wave functions \ips) = |0000) and \iPt) — 
|000, -1), we obtain 
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FIG. 5: (Color online) The form-factor J^sT{q\\) for \/a% 
0, 1, 3, as well as 10^ J^sT{q\\) for A/as = 10^. 



where 7^ and Am, with m = 0,1, are the variational 
parameters found in Appendix [X] Note that J-'st is in- 
dependent of the in-plane angle of the emitted phonon, 
due to the the circular symmetry of the confining po- 
tential. In the absence of Coulomb interaction (7™ = 
and A = A), Eq. (|98p reduces to the following expression 
Wa*B = 0) 



^ST {q\\] 



8m*uj 



(99) 



where we used that A — y^2h/m*LLj and A = y^h/2'm*uj. 
Note that the function J^st {q\\) in Eq. ([99| has a single 
scale, g|| ~ ^jAmFuj/h. For strong Coulomb interaction, 
however, we find that Tst {q\\) has two scales; an addi- 
tional scale arises due to the Coulomb interaction and 
it is given by g|| ^ '^/{r), where (r) = A(A/2ag)^/^. In 
the limit X/a*g > 1, we find Tst {q\\) from Eq. ^ 
for both g|| ^ ^/{r) and g|| ^ 1 /('"), and then match 
the two asymptotes into the following crossover function 

{\/a*g » 1) 



^ST {q\\ 



(1\\0-B 



Jl 



(l + v/3)fi,2 



(100) 



where (r) = X{X/2a*sy/^ and Ji(x) is a Bessel function of 
the first kind. We checked numerically that Eq. (jlOOp co- 
incides with Eq. ((98)) in the whole range of g|| for A/a|j ^ 
1. By comparing Eq. (fTUO)) to Eq. ([M]), we find that the 
maximal value of Tst {q\\) scales like oc (a|j/A)^/® for 
strong Coulomb interaction. In Fig. [51 we plot J-'st {q\\) 
as a function of g|| A for X/a*g — 0, 1, 3, which shows that 
J- ST {q\\) gets suppressed by the Coulomb interaction for 
A/a|j > 1. In addition, we also plot lO^JFgr (gy) for 
X/a*^ — 10^, which shows that J^st {q\\) oscillates on the 
scale of Ag|| = l/(r), see also Eq. (jlOOl) . 

Next, we find the relaxation rates Fq,^ by substituting 
Eq. (pS)) into Eq. (PTjl and averaging over phonons at 
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thermal equilibrium. We obtain 



2 2 

i Q/3 = a^Qp 



^ TA 



(w/3a) 



where = {Efj —Ea)/fi is the transition frequency, the 
coefBcients are given in Eq. (|96p . and 



(101) 



ff^(u;) = (1 + iV^) 

r£f(^.) 



nhpcsl 



tt/2 



M 



TsT ( ^-sini^j f ^cosi9 



'14' 




ir/2 



TsT ( — sini? j F f — cost? 



sint9, 

2 



1/7 \ f W 

J^ST ( — sin I? — cos "d 

S2 / VS2 



sin^ tJcos^ 

(8 cos'' t? + sin'' d) sin^ ?9, 



(102) 
(103) 
(104) 



where iV^ = {e^-^/T _ The subscripts LA and 

on the left-hand side in Eqs. (|102p - (ll04p correspond to 
the dot electrons interacting with either a longitudinal 
acoustic {LA) or a transverse acoustic ( TA) phonon. The 
superscripts give the mechanism of the electron-phonon 
interaction: deformation potential (-DP) or piezoelectric 
{PE). In Eq. (|101[) . we have taken into account only the 
acoustic phonons. The optical phonons can be neglected 
due to their large excitation energy [Hlulo ~ 37meV and 
fvjjTO ~ 34meV, at T = in GaAs). For the acoustic 
phonons, we used a linear dispersion relation, — sjq, 
which requires that the transition energy hujpa be much 
smaller than the Debye energy (« 30meV). In prac- 
tice, the linear dispersion relation is fairly accurate for 
hujLA < 23meV and Hujta < 8meV, in GaAs. In this 
material, the speed of sound is si = 4.73 x 10^ cm/s for 
LA phonons, and S2 — S3 — 3.35 x 10^ cni/s for TA 
phonons. As mentioned above, we take into account both 
the deformation-potential (Sq = 6.7 eV) and piezoelec- 
tric (/ii4 = — 0.16C/m^) mechanisms of electron-phonon 
interaction. Finally, we assume a dielectric constant 
K = 13.1 and a crystal mass density pc = 5.3g/cm'^. 

We consider further the relaxation rates: r_|__Oj Tq,-, 
and r_(_^_. The other relaxation rates are either zero, in 
this approximation, or can be found from the detailed- 
balance principle in Eq. (|92|) . Note that the states |'I'±) 



and l^'o) are linear combinations of the singlet \'^s) and 
triplet |5't±) states, as given in Eq. ([M)) . Also recall 
that the triplet state |\E'to) decouples for these approxi- 
mations, and hence it is long lived. To analyze the re- 
laxation rates F+ Oi Tq,-, r-|_^_, we first consider the fac- 
tor a^a| in Eq. (|10ip separately. As before, we assume 
A± <^ Ez, or equivalently (r) <C A50. In addition, we 
also assume that A+ and A_ do not differ much from 
each other, i.e. A± ^ A^/i?^- Under these assumptions 
and for Ez = —\Ez\^ we obtain the following expressions 
for alal, 

2ao4 = j2 1 4A2 1 1 ± / , ' (105) 

4a^«2_ ^L^ J+ j L J \ (log) 

where J± = Et± — Es (or J± — Ets ± L^z)- In case of 
Ez — \Ez\, one should replace J± — > and A± 
in Eqs. (jlOSp and (|106p . In Fig. ^ we plot a^a^ versus 
B across the singlet-triplet crossing in a model specified 
in the caption of the figure. 





Ets > \Ez\, 
Ets ~\Ez\-> A+ 


Ets ~ \Ez\, 
\Ets ~ \Ez\\ < A+ 


\Ez\ > Ets > -\Ez\, 
\Ez\ - \Ets\ > A± 


Ets ~ -~\Ez\, 
\Ets + \Ez\\ < A- 


Ets < -\Ez\, 
-Ets - I-BzI > A_ 


2 2 

aoa+ 


AI/{Ets - \Ez\f 


1/4 


AI/{Ets - |£z|)' 


< {A+/2Ezf 


AlAl/{Els-El)' 


2 2 


^lA^^KEls - El)' 


< {A-/2Ezf 


AI/{Ets + |Sz|)2 


1/4 


AI/{Ets + \Ez\)' 


2 2 


/\'-/{Ets + \Ez\)'' 


< {A-/2Ezf 


AlAi/iE^s-El)' 


< iA+/2Ezf 


AI/{Ets - \Ez\)' 



TABLE I: The factor a^a^ in different regimes of Ets- We assume here Ez = — |-Bz|, however, the table can also be used for 
Ez = \Ez\ after replacing A± A^. 
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FIG. 6: (Color online) (a) The factors aoa+ (solid line) and 
aoi- (dashed line) in Eq. pOSp as functions of the magnetic 
field Bz = B. (b) Similarly, the factor aiat in Eq. ((106)) as 
a function of B. For the plot we used the following model 
(arbitrary units); A+ = O.lSi^zl, A_ = Q.l\Ez\, Es = 0, 
-Et± = IQ ~ 2B ± Ez, and Ez = -0.2B. The exact solu- 
tions for a^a^ (dotted line) coincide with the approximated 
ones (see Eqs. (|lU5p . (|lU6p ) on the scale of the plot. Note the 
different scales on the left and right panels. 



From Eq. (fT05)) . we find that 



has a Lorentzian 
line-shape as a function of Ets in the vicinity of Ets — 
ztjiJ^I, see also Fig. IUJa). The maximal value of a^ai^ 
approaches 1/4 in the vicinity of Ets = ^\Ez\- Around 
this point {Ets = the spin and orbital degrees of 

freedom of the two electrons are strongly intermixed be- 
tween themselves. As a consequence, the relaxation rate 
r+.o (ro,-) is independent of the spin-orbit interaction 
strength, and it is given by (l/4)r(2A±/ft) according to 
Eq. poip . The quantity T{'w) depends on the type of 
bath that is coupled to the electron charge. Whether an 
appreciable relaxation takes place or not depends thus 
on the value of T{w) in the limit w 0. As we show 
below, r(w) due to the electron-phonon interaction is 
generally suppressed for w < s/X for all mechanisms of 
electron-phonon interaction. Therefore, the resonant en- 
hancement of spin relaxation at Ets = ^\Ez\ can occur 
due to phonons only in structures with A± ~ hs/X. How- 
ever, A± ^ hs/X in typical GaAs quantum dots used 
for studying spin relaxation, and therefore the phonon- 
induced spin relaxation is suppressed at Ets = ^\Ez\ in 
these structures. Finally, we note that a^at is strongly 
suppressed for Ets > \Ez\, and a§a^ is strongly sup- 
pressed for Ets < —\Ez\- In both cases a^a^. is sup- 
pressed owing to the selection rule of Hso^ which al- 
lows only a single spin to be flipped in a first order pro- 
cess. Thus, for example, the state l^"-) at Ets > \Ez\ 
consists mostly of \^tJ) and has an extremely small 
(second order in Hso) admixture of |^'t+)- Similarly, 
|*o) ~ |*T+) +C'(-ff|o)l*r_)- As a result, a scalar fluc- 
tuation, like C/ph(i), does not connect |^-) and I^Po) in 
the flrst order of Hso ■ 

Next, we consider the factor a\a^_. From Eq. (|106p . 

we find that a^a?_ = A?_/Ji on the left-hand side 

~ r2 



side {Ets < —\Ez\) of the singlet-triplet crossing. a^a?_ 
as a function of Ets has two maxima, at Ets ~ ^Ez, 
see Fig. Eljb). The upper bounds on these maxima read: 
a^a^ < A^/4£'| for the maximum at Ets ~ \Ez\, and 
a\a^ < A^/4i?| for the maximum at Ets ~ — \Ez\- In 
the region between the maxima (— jE'^l < Ets < \Ez\), 
the factor a^a^ is strongly suppressed owing to the se- 
lection rule of Hso, as discussed above. We summarize 
our results for a'^a'^ and a^a^ in Table HI for a number 
of limiting cases. Except for the two regions of avoided 

2 „2 



/5 



crossing, at Ets = ^Ez, one of the three factors a^a 
is always suppressed compared to the other two. As a 
result, for most values of Ets, there are only two relax- 
ation rates of interest, which correspond to transitions 
between the triplet states |^t±) and the singlet l^s). 

We turn now to analyzing the terms in the square 
brackets in Eq. (I101|) . The spin-orbit interaction en- 
ters here indirectly, via the transition frequency Lopa- 
The latter is considerably influenced (renormalized) by 
spin-orbit coupling only at the points of avoided cross- 
ing {Ets = ±Ez). Each rate t{w) in Eqs. (fT02l) - (fT04l) 
is thus independent of the spin-orbit coupling and can 
be studied separately as a function of the transition fre- 
quency w. For simplicity, we assume T = 0, which allows 
only spontaneous emission of phonons to occur {N^ = 0). 
In the low frequency limit, w <^ ^/(r), wc obtain from 

Eqs. (nni-dMi) 



(107) 



Such a frequency scaling is valid for all quantum dots that 
have conflning potentials with in-plane inversion symme- 
try, U{x,y) = U{—x,—y). For such quantum dots, the 
form-factor TsT{q\\) is real and scales as TsT{q\\) oc at 
q\\ l/{r). In contrast, for asymmetric quantum dots, 

the terms ^^^(w) oc and i'i,A{TA)('^) P°^" 
sible in the leading order of w. For our quantum dot, 
with U{x,y) = m*LUQ{x'^ + j/^)/2, we obtain explicitly 
{w < s/(r)) 



T^taM 



, A(A/a|), 
?^^^^-I^±!^Ao{X/a*s) 



(108) 
(109) 
(110) 



where Ao{X/a*g) < 1 is a suppression factor, due to the 
Coulomb interaction, given by 



An{X/al 



-^1,0 



A' 



n+4 



(l + v/l + 7i)A?-(l + ^/^)A2j , 

(111) 

with 7m and Am, for m — 0, 1, being the variational pa- 
rameters found in AppendixIXl and Ai^o given in Eq. 
{Ets > \Ez\) and a^a^ — A^/J^ on the right-hand Note that An{X/ag) depends only on the ratio X/a*g. 
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FIG. 7: (Color online) The factors Aoi^/a's) (solid curve), 
A2{^/a*g) (dashed curve), and A5{X/a*g) (dot-dashed curve) 
as functions of the Coulomb interaction strength X/ag. Inset: 
the factor A5{X/a*g) on a larger scale. 



We plot Ao{X/a*^) as a function of X/a'g in Fig. [7] (solid 
curve). In limiting cases, we obtain 



Ao{x) 



(112) 



Thus, the Coulomb interaction has a strong effect on the 
spin relaxation rates. 

Assuming Ez ^ Ets ^ hs/{r) and Ets,—Ez > 0, 
we find the triplet-to-singlet relaxation rates (rgT-i- = 
r+o and TsT- = r+_) 

IZire'^h'i.E'^aElX^ . . , . . 



where we introduced two effective spin-orbit lengths. 




(114) 



In Eq. (|113[) . we retained only the piezoelectric inter- 
action with transverse phonons (F^^), because it gives 
the largest rate for GaAs quantum dots. The suppres- 
sion factor A2 (A/a|j) is shown in Fig. [7] (dashed curve). 
For weak Coulomb interaction, we have A2 (A/a|j) ~ 
1 — 0.203 (A/a^) , and, for strong Coulomb interaction, 

we have A2 (A/a|j) = (4/9) (A/2a|j)~^/^ Thus, in this 
case, the effect of Coulomb interaction is to reduce the 
singlet-triplet relaxation rates. This reduction is despite 
the fact that the Coulomb interaction enhances the spin- 
orbit effects in the quantum dot, e.5., the splittings A± in 
Eq. (|65p . see also Fig. [1] The enhancement of spin-orbit 
effects occurs due to an effectively larger quantum-dot 
size, Ai,o ~ {r). Note that Ai^q/A = 1 + (V^/8)(A/a|j), 
for \/a*B < 1, and Ai,o/A = (A/2a|j)i/3^ for X/a*^ > 1. 



therefore (Ai^o/A)^ can be viewed as an enhancement 
factor and ^0 (A/a|j) as a suppression factor. The sup- 
pression factor ^0 (A/a|j) originates from the form-factor 
!Fst{<1\\) in Eq. ([55]) . which is sensitive to the difference 
|V's(ri,r2)|2 - |V'T(ri,r2)|2, see Eq. 1^. The Coulomb 
interaction makes the singlet and triplet charge profiles 
look more similar to each other, and thus, reduces the 
relaxation rates. An additional reduction of relaxation 
rates stems from ETS^ which enters Eq. (|113[) . The de- 
pendence of Ets on the Coulomb interaction strength is 
studied in Section Hill sec Eqs. ^ and (jB2|l . 

Assuming now \Ets—\Ez\\ ^ A+ and \Ez\ <C 
hs/{r), we obtain three relaxation rates, 



_ 13V2nff4^\EzfX^ 



, ^ 26TTe^hl^\Ez\^X'' 



^2 (A/a 



(115) 
(116) 



where the factor ^5(A/a|j) > 1 is shown in Fig. [7] 
(dot-dashed curve). Note that ^5 (A/a^) = 1-1- 
{3y/^/8){X/a*g), for X/a*^ < 1, and A (A/a|j) = 
(4/9) (A/2a|j)^/^ for X/a% > 1. The relaxation rate 
F+.o in Eq. (|115p is enhanced by the Coulomb interac- 
tion, via the factor ^5(A/a|j) > 1, whereas the relaxation 
rates Fq,- and F4.^_ are suppressed, due to the factor 
A(A/a|j) < 1, like in the previous regime. The enhance- 
ment due to the factor A5{X/a*g) is not significant for a 
realistic structure, however, it confirms our picture that 
there are two competing effects of the Coulomb interac- 
tion on the singlet-triplet relaxation, as described above. 
In the case of A(A/a|j), the increase of the quantum 
dot size ((r) « Ai 0) wins over the suppression due to 

/o'oW-/o'iW- 

Considering further Ets ^ Ez ^ hs/ (r), wc obtain 



_ 13ne^hl,\EzfX^ 

^ +^0 - ioogn..2fcf;^ 7A2 -^2(A/ag) 



13860K2ft6pcS2A^ 



13Tre^hl^\Ez\^X^ 
13860 K^h'^PcslA'i 



A(A/a| 



(117) 
(118) 



Also note that A2 {X/a*g) 



(Ai,o/A)Mo (A/a|j), and 



where ^2(A/a|j) was analyzed above. Again, we see that 
the Coulomb interaction reduces the relaxation rates in 
this regime. 

On the triplet side of the crossing {Ets < 0), one can 
use Eqs. (fTT5)) and plB)) with F+^o ^ Eq _ and A± 
for the regime Ets « -\Ez\ < fis/(r), and Eqs. (fTT^ 
with Ts.T± ^T±,s and Ets ~^ \Ets\ for the regime 
Ez < Ets < fis/{r). 

So far, we have considered only the regimes of small 
transition frequencies, w <^ s/{r). The quantities F(?«) 
in Eqs. p02p — (|104p are functions of the transition fre- 
quency w, and contain several scales of w: s/(r), s/X and 
s/d. In principle, s/(r) can be much smaller than s/A, 
provided the Coulomb interaction is strong enough. In 
the regime s/(r) ^ w <C s/A, we find oscillations of F(i/7) 
as function of w with a period Aw ~ s/{r), see Fig. |51 
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FIG. 8: (Color online) Schematic of I'LAi'w), showing the 
relevant scales oi w: s/X, s/d, and s/{r). The curve was 
obtained using the form- factor in Eq. (|100|) . for a fairly large 
value of the Coulomb interaction strength, X/ag = 10'^. 
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FIG. 9: (Color online) The factors Bi, B2, and B3 as functions 
of the Coulomb interaction strength X/a*g. Inset; the factors 
in the main figure multiplied by (Ai.o/A)^. 



These oscillations stem from the oscillations of TsTiqw) 
in Fig. [51 Note, however, that a fairly large strength of 
the Coulomb interaction is required to observe these os- 
cillations. The reason for this is that (r) = X{X/2a'^y^^ 
scales weakly (power 1/3) with A/a|j ^ 1. In a typical 
GaAs structure with X/a*g ~ 5, there is not sufficient 
room between the scales s/{r) and s/A for these oscilla- 
tions to occur. 

In the regime s/X ^ w ^ s/d, we find 



r 



2a (w) oc w, 

r^f (w) cx w 



-3 



(119) 
(120) 
(121) 



These scaling laws are universal in the sense that the 
power of w does not depend on the dot confining po- 
tential, provided A is understood as the smallest lateral 
scale. The phonon of frequency w is emitted at an angle 
■d ~ arcsin(s/Aw) to the z-axis. From Eqs. (|102p - (|104|) . 
we obtain 



2nhpcSi\'^ 



Bi{X/a*s), 



17287re2ft.2^s 



K?hpcX^W^ 



-B2{X/al), 



K^hpcX'^W^ 



BsiX/a* 



(122) 
(123) 
(124) 



where the suppression factors S„(A/a|j) are shown in 
Fig.H We note that Bn ^ {X/a*^)-'^/^ at X/a*g > 1. 

Next one can consider again different regimes given in 
Table [T] and using Eq. (jlOip obtain relaxation rates for 
different mechanisms of electron-phonon interaction. We 
skip this part here and consider only the experimentally 
most relevant regime Exs ^ Ez,hs/X. For GaAs, we 
estimate hs2/X « 22 peV for A = 100 nm. Retaining only 
the deformation potential mechanism, we obtain using 



Eq. dual, 



772 ri2 \2 



IGTrfi^PcslETs^lX' 



:Si{X/a*s), 



(125) 



for Ets > 0. On the other side of the crossing {Ets < 0), 
we have the same expression, but with Ts,t± ^ ^t±.s 
and Ets \Ets\- Finally, we note that {Xifi/X)'^Bn 
is a monotonically decaying function of A/a|j, see inset 
of Fig. [9l and thus the Coulomb interaction reduces the 
relaxation rates. 

In the regime w S> s/d, the scaling of T(w) with w is 
non-universal, in the sense that it depends on the quan- 
tum dot lateral and transverse confinement potentials. 
Considering the quantity 



max 

1? 



S2 



— sm v \ }' [ — cos v 



S2 



(126) 



one can find a line ??(?«) in the {w, z?)-plane such that 
for a given w the product of form-factors is maximal at 
i9 — ^{w) . Then, except for the special cases: ■d{w) — 
0, 7r/2, the phonon emission is peaked at 'd{w). The scal- 
ing of r(w) with w follows the square of the expression 
in Eq. p26p . up to some prefactors which are power- 
law in w. In the case of our harmonic lateral confine- 
ment, the form- factor J-'sTiQ\\) decays Gaussian-like at 
g|| ^ 1/A. Provided the transverse form- factor ^(q^) has 
not a stronger decay law, then the least-suppression line 
is i){w) = 0. The phonons are then emitted at an angle 
close to 0,7r (more precisely at ~ arcsin (s/Aw)). The 
scaling of T{w) with w follows the scaling of |F(w/s)p, 
up to prefactors which are power-law in w. In the other 
extreme case, when F{qz) has a stronger decay law (at 
Qz » l/d) than J^st{q\\), then 't}{w) ^ arccos(s/dw) and 
T{w) (X |JFst(w/s)P. Thus, the regime of large transi- 
tion frequencies (w 3> s/d) allows one to learn about the 
dot confinement, however the singlet-triplet relaxation 
rates are strongly suppressed in this regime. For GaAs, 
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FIG. 10: (Color online) Relaxation rates F+^o, r+,-) a^nd Fo,_ 
as functions of an in-plane magnetic field B — (0, -By , 0). An 
avoided crossing of the singlet 15) and triplet |r+) occurs at 
Ets = -Ez « 0.75 meV (By « 8.2 T). We used parameters 
typical for a GaAs quantum dot: hujo = 2meV, X± = 8 /im, 
d = 5nm, g — —0.44, and k — 13. Inset: Schematic of level 
crossing for applying an in-plane magnetic field. 



we estimate that the singlet-triplet energy Ets needs to 
be larger than hs2/d « 0.4 meV at d = 5nm, for this 
regime to be accessible. 

To summarize this section, we present the calculated 
relaxation rates as functions of the magnetic field B. For 
an in-plane magnetic field B ~ (0,i?||,0), we show the 
relaxation rates in Fig.fTUl Only a single avoided crossing 
(at Ets = —Ez) is present in this case, see inset of 
Fig. [TOl The states and ^5 are strongly intermixed 
with each other, which results in a large relaxation rate 
between these states. For the parameters used in Fig.lTUl 
the relaxation between ^fT^ and 'i's (see F+.o in Fig. fTU]) 
is about 10^ times stronger than the relaxation between 
^'t_ and ^'5. At Ets — ~Ez, the relaxation rate F+^o 
has a dip, because the phonon-emission energy hw = 
2A+ is very small (A+ ^ hs/X) at this point. We note 
that the value of r+_o in the dip depends strongly on A+, 
as it can be seen after rewriting Eq. (jllSp as follows, 



+,0 



(127) 



The sensitivity of A+ to the magnetic field direction, 
illustrated in Figs. [H and [21 can be used here to tune 
F+ by orders of magnitude. It is also important to note 
that at this point {Ets = —Ez) the eigenstates I^I'o) and 
1 4'+) are approximately given by 



l*T+)) 
l*T+)) 



(128) 
(129) 



where 1*5) = \S)\tlJs) and \'^t+) = \T+)\^pT), with 
\i^S{T)) being the singlet (triplet) orbital wave func- 
tion. In Eqs. (|128p and p29p . we neglected terms cx 



0(A/Aso)|^T_)- Thus, at this degeneracy point (and in 
its vicinity), there is no spin selectivity in coupling to 
potential fluctuations, because the spin is not a quan- 
tum number here, see Eqs. (|128p and (|129p . While the 
phonon-induced relaxation is suppressed at small fre- 
quencies, the coupling to some other degrees of freedom 
in the solid state {e.g. l//-noise, noisy gates, i?-field fluc- 
tuations, particle-hole excitations on the Fermi sea, etc.) 
can dominate the relaxation between \ and |vE'+). The 
relaxation rate is given by the orbital (charge) relaxation 
rate at the energy hw = 2A+. Furthermore, since the 
states (|128p and (|129p arc difhcult to discriminate and 
prepare by the currently used techniquesj^iii^i^ it is 
more convenient to study relaxation at the avoided cross- 
ing point by means of the Ramsey fringes. Then the re- 
laxation rate in Eq. (|127p appears as a decoherence rate. 
The required rise (fall) time of pulses for tuning to the 
degeneracy point should be smaller than ttH/ A_|_ « 200 ps 
for A = 100 fieV, which is feasible experimentally (note 
that Ets can be tuned by electrical gates on sort time- 
scales, while keeping Ez constant). 

Far away from the dip in Fig. [TOl say to the left {\Ez \ < 
Ets) 7 the relaxation rates read 



27rh^p,slXHETsT\Ez\)' 



(130) 



where we assumed the regime s/A ^ w <C s/d, with 
w = {Ets T \Ez\)/fi, and took into account only the de- 
formation potential mechanism. The latter approxima- 
tion is valid only for large transition frequencies w, typ- 
ically > 107re|/ii4|s^/KSoA. At smaller w the piezo- 
electric TA mechanism gives the largest contribution to 
the relaxation rates. Thus, the rate Tqt-^- turns to 



96TTh^e^hl^s2AlB3 
>i^pAHETS-\Ez\)' 



(131) 



as Ets ~ \Ez\ is decreased, but still remaining in the 
regime s/X < {Ets ~ \Ez\)l^ < sjd. 

For a magnetic field perpendicular to the 2DEG plane, 
B — (0,0, i?j_), the relaxation rates have a richer be- 
havior, because more than one avoided crossing can oc- 
cur in this case, see for example Fig. [1] Furthermore, 
A_|_ and A_ can strongly differ from each other, due 
to their different dependence on the i3-field direction, 
see Eq. (|65|). In Fig. [11] we plot the relaxation rates 
as functions of the magnetic field which is tuned 
across the singlet-triplet crossing. Two avoided crossing 
occur at B w 0.85 T {Ets = \Ez\) and at B w 1.1 T 
{Ets — ~\Ez\)- We choose A+ — 2A_, which means 
that A_ = 3A+, see Eq. (jHS]) . The relaxation rates in- 
crease strongly as one approaches the neighborhood of 
the singlet-triplet transition from either side (maxima of 
rates correspond to the optimal phonon emission condi- 
tion: transition energy on the order of hs/X). In con- 
trast to the case of in-plane field (see Fig. [TO)) , both re- 
laxation rates are strongly affected here. Although the 
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FIG. 11: (Color online) Relaxation rates F+^o, r+,-j and Fo,_ 
as functions of the transverse magnetic field B = (0,0, -Bx). 
Two avoided crossings occur at Ets = ^Ez, as shown in 
Fig. [T] We use the same parameters as in Fig. IIUI except 
for A+ = A-/2 = 8 /.im. For these parameters, we obtain 
Ets ~ 0.75 meV at = 0, and Ets = at Bx ~ 0-95 T. 



overall behavior of the rates in Fig. [TT] is more compli- 
cated than in Fig. [TOl the behavior of the rates in the 
vicinity of an avoided crossing is not much different. For 
definiteness, let us consider the avoided crossing on the 
left-hand side, occurring at B w 0.85 T or equivalently 
at Ets = —Ez = \Ez\- There, the rate F+^o exhibits a 
dip as a function of S^, whereas the rate r-|_,_ (Fo,_) de- 
creases (increases) abruptly at the point of avoided cross- 
ing. The occurrence of the dip is explained, as before, 
by the suppression of phonon emission at small energies 
hw = 2A+ ^ hs/\, see also Eq. (|107l) . This behavior is 
similar to the behavior of the rate F+.o in Fig. [TOl The 
abrupt dependence of the rates F_|_^_ and Fo,_ on B±^ at 
the avoided crossing is due to the fact that the transition 
occurs between a distant level l^"-) « \^t^) onto two 
closely spaced levels and I'I'o), which change from 

l^-^-) « l^-s) and l^-o) ~ |*T+) on the left-hand side of 
the avoided crossing to \'^+) « |*t+) and |^'o) ~ \^ s) 
on the right-hand side. Obviously, the transition between 
|\1/t_) and \^t+) is forbidden due to the spin selection 
rules, and therefore, the rate Fo,_ (F+__) is suppressed on 
the left-hand side (right-hand side) of the avoided cross- 
ing. As one sweeps the magnetic field across the avoided 
crossing the rate F+,_ (ro,-) decreases (increases) in 



the neighborhood — A_|_ < Ets + Ez < A+. This in- 
crease/decrease appears abrupt on the scale of the plot, 
because A+ is very small for GaAs (A± ^ Ez)- Note 
that the same abrupt-change behavior is present for the 
rates F+^_, Fq,^ in Fig. [TOl The second avoided crossing 
(at B m 1.1 T) in Fig. [TT] shows an analogous behavior 
with the rate Fq,- exhibiting a dip and the rates F+^o 
and F_|_._ interchanging their roles. We thus conclude 
that the physics of the singlet-triplet relaxation for per- 
pendicular and parallel B-fields is very much alike, al- 
though the case of perpendicular field is richer because 
two avoided crossings take place. As for tilted magnetic 
fields, we did not find any qualitative difference as to 
the case of a perpendicular field, as long as two avoided 
crossings were present in the considered S-field range. 



B. Triplet-triplet relaxation rates 

Away from the singlet-triplet crossing, the triplet con- 
sists of three levels symmetrically split by the Zeeman 
energy, Et± = Et ± Ez and Etq = Et- The transi- 
tion between |r_) and \T+) is forbidden at the first or- 
der of Hso, as we mentioned already in Eq. ([QS]) . The 
transition between |To) and \T±) is allowed and occurs 
due to virtual processes with going to states outside the 
singlet-triplet subspace (see Fig. |4l[a) and Sec. IV C|) . At 
the singlet-triplet (avoided) crossing the states \T±) mix 
with the singlet state (see Sec. IV Bl) . We are, therefore, 
entitled to consider the relaxation of |To) onto the eigen- 
states 

|*a) =aa|*s)+fca|*Tj+Ca|^'T_), (132) 

where a = 0, ± labels the eigenstates and the coefficients 
aa, ba, and Cq were found in Sec. lVBl Next, we use Ui^t 
given in Eq. with given in Eq. In Eq. ((75|l . 

we retain only the first term (see Fig. 3] and discussion 
in Sec. IV Cp . We evaluate first the matrix elements of 
U\at with respect to the unmixed states. As a result, we 
obtain 

(^T±|f/i„t(t)|*To) = T^Ez{X T lY) ■ n{t), (133) 

where f2 = {^x', ^y' , 0) is a dimensionless quantum fluc- 
tuating field. Its components are given by 



^fe'Afq,.fT(g||)F(g,) + 



(134) 



where q = (q||,'7z) and q|| = {qx'jQy'). In Eq. (|134p . !FT{q\\) denotes the following form-factor, 

^T(q||) = (^T|e"iii^cos(q||r/2)|V'T) 



(135) 
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We evaluate TT{q\\) for the triplet state |i/'t) = |000, — 1) 
and obtain 

(136) 

where 71 and Ai are the variational parameters for m = 1 
studied in Appendix |^ Note that Tt is independent 
of the in-plane angle of the emitted phonon, due to the 
circular symmetry of the confining potential. 

Next we use the states in Eq. (|132p and the matrix 
elements (I133P to evaluate the relaxation rates according 
to Eq. (HI]), 

TaTo - (1 - al)T^X + 2Re[6„c:]r«^ + 21m[6„c:]r(f^^ . 

(137) 



The first term in Eq. (|137p gives the triplet-triplet re- 
laxation without accounting for interference between the 
two processes in Eq. (|133p . The other two terms in 
Eq. (|137p are due to the interference, and are therefore 
relevant only close to the singlet-triplet crossing. Away 
from the singlet-triplet crossing the first term alone gives 
the triplet-triplet relaxation rates. To shorten notations, 
we consider further F^^"^, where either a = Tg or /3 = Tq. 
We find 



Q/3 



E 



r.Tr/2 



sin-d F 



cos-d 



(138) 



where topa = {E/3 — Ea)/h is the transition frequency. 
The effective spin-orbit length A50 is defined as follows 



'1 



Aso 



A' 



(139) 



In Eq. (jl38p . the summation runs over the acoustic 
phonon branches, with j = 1,2,3 corresponding, respec- 
tively, to longitudinal and two transverse polarizations. 
The quantities and in Eq. (|138p are introduced as 
follows 



rate. They are, therefore, not positively defined. We 



find for r^^^^ and T^l^^ similar expressions as for Ej^^^^ in 
Eq. (|138p . but with a different geometrical prefactor. To 
obtain E^^'^^ and E^^^^^ from Eq. (I138p , one has to replace 
the factor 1/A|q by 
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Al 



Iz sin 2if' , 



tF. (142) 



(143) 




dcj> cos" 



(140) 



After integrating over the phonon angle we obtain ex- 
plicitly 



^0, S2 — ^3 — 0, 
3-\/27r/ii4 



K 

V2T:h 
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■ sin d cos i?. 



sin2i?, 



-(3cos2^?- l)sini?. (141) 



We note that the explicit form of /32,i? and Ps^^ depends 
on the choice of transverse phonon polarizations, whereas 
i9 + 1? is independent of this choice. 

The quantities F^"';^^ and F^^^^ in Eq. p37p do not have 
the meaning of rates, but rather of corrections to the 



respectively. Here, ip' is the azimuthal angle of I. 

Next, we evaluate the product 6q,c* using Eqs. ([71]) 
and (|63p . The resulting expression (with separated real 
and imaginary parts) reads 



\2 „2 
8£'qT+ EaT^ 



A 



so 



(144) 



where EaT± = E^ — Et^ . Using Eq. (|144p and the pro- 
portionality of F^ 
rate in Eq. (flSTP 



portionality of F^"'^^ and F^^^^^ to Y^^Tq ' obtain for the 
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(0) 
aTo 



= 1 



'^EaT. EaT- 



A'"* 
^^so 



A"4 
^^SO 



(145) 

The latter equation allows one to evaluate the relaxation 
rate ToTo knowing the "bare" rate F^"^^ , which in its turn 
can be calculated from Eq. (|138p . The physical meaning 
of F^°^^ becomes clear if one considers taking the avoided 
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crossing splittings A± to zero, while keeping the virtual 
transitions {weak matrix elements) finite. Then, there 
exist only two rates, Tt_To and ^ToT+- The reverse rates 
are trivially given through the detailed balance equation 
([92)1 . Due to spin-symmetry, these rates are equal to 
each other, Tt_To — ^TqT+, and can be considered as 
a function of the transition frequency. This function is 
nothing but the rate T^^^ in Eq. (|138p ; note that the only 
dependence on the indices a and f3 in t'"^^ is through 

the transition frequency w^a- Therefore, T^^^ has the 
meaning of the "bare" triplet-triplet relaxation rate. 

The "bare" rate r^^"^ is closely related to the one- 
electron spin- flip ratci^i^ To illustrate this, we briefly 
outline the derivation of the one-electron rate. It is 
straightforward to verify that, for one-electron case, 



Eq. (|133p takes the following form 

(*TlC^int(OI*i> = -iEziX - iY) ■ n^it), (146) 

where fli{t) is the field acting on a single electron. An 
expression for fli{t) can be obtained from Eq. (jl34p by 
replacing JFt(<Z||) ^ ^•^o(9||)- Here, J-'o{q\\ ) is the single- 
electron form-factor, given by 

- (^ole^'^ii'^M^o) = e-^«^^ (147) 

where we assumed that the dot is in its ground state 
Mri) = A-V-i/2exp {-rl/2\l), with Ad - y/h/m*uj. 
The spin-flip rate is then given by the following expres- 
sion (assuming Ez = — \Ez\ and denoting loz = \Ez\/fi) 



- 



dd — ^ 



To { — sini? ) F ( — cosi9 



(148) 



Note that Eqs. (|138p and (|148p are very similar to each 
other, differing only by \Tt {q\\) |^ ^ \ (^ll) 1^ gO" 
ing from Eq. to pig]) . The form-factors Tt {q\\) 

and J-Q ((Jii) differ significantly only at large strengths of 
the Coulomb interaction. In the absence of the Coulomb 
interaction, Eq. (|136p gives 
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(149) 



which differs from Eq. (|147p only by the prefactor in the 
parentheses. Note that A = \/2\d and ii^i(— l,l;x) ~ 
1 — X. For strong Coulomb interaction, we find the fol- 
lowing (approximate) crossover expression for JFt('?||), 
by matching the limits of small and large q|| , 



Tt {q\\) = Jo 



<l\\{r) 



(150) 



Equation differs from Eq. pi7)) both by the pref- 

actor and the argument of the exponent. We note that, 
since the form- factor Tt {q\\) is integrated over its ar- 
gument in Eq. (|138p . the resulting rate is, by order of 
magnitude, equal to the one-electron spin-flip rate in 
Eq. p48p . except for the case of strong Coulomb inter- 
action. In the case of strong Coulomb interaction, a new 
frequency scale, uj ^ s/{r), emerges in the triplet-triplet 
relaxation, analogously to the singlet-triplet relaxation 
studied in Sec. IVI Al 

Finally, we note that for weak spin-orbit coupling the 
interference terms in Eq. (|137p can be neglected, since 
they are smaller by at least one power of X/Xso as com- 
pared to the first term. In this case, the triplet-triplet 



relaxation is governed by the transitions r_ ^ Tq and 
To <-!■ r+, both of which occurring at the same rate. 
For weak spin-orbit coupling, it is also hard to discrimi- 
nate between eigenstates and bare states at the points of 
avoided crossings in usual spin relaxation experiments. 
Therefore, for the accessible regions, the factor 1 — is 
either zero or unity. Thus, we arrive at a single non-zero 
rate ^t^To = ^ToT+ , which is obtained by evaluating T^^^ 
in Eq. (|138p at the transition frequency = \Ez\/h, 
for Ez — —\Ez\- The reverse rate Tt^To = ^TqT^ is 
obtained by multiplying the result by exp {—\Ez\/T). 



VII. CONCLUSIONS 



In this paper, we studied spin relaxation in a quan- 
tum dot with two electrons. The spin relaxation occurs, 
as in the case of single-electron quantum dots, due to a 
combined effect of the spin-orbit, Zeeman, and electron- 
phonon interactions. We focused on the vicinity of a 
singlet-triplet crossing, which occurs with applying an 
orbital magnetic field, considering thus relaxation be- 
tween states of the singlet-triplet subspace. We found 
that the spin-orbit interaction induces avoided crossings 
between the singlet \S) and the triplet states \T±) as the 
quantum dot is tuned across the singlet-triplet crossing. 
At the same time, the degeneracy at the crossing of the 
singlet \S) and the triplet state |To) is not lifted at the 
first order of the spin-orbit interaction. This result has 
the implication that the spin relaxation occurs efficiently 
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only between the following states: 

S^T+, S^T^, (151) 
To ^ r+, To ^ r_. (152) 

Close to an avoided crossing point, the singlet and triplet 
{T±) states are superimposed with each other and form, 
thus, states of no definite spin. We find, nevertheless, 
that the phonon emission is inefficient in this case for typ- 
ical GaAs quantum dots, because the transition energy 
is too small for the spin-phonon coupling to be relevant. 
As usual, the strongest coupling to the phonons occurs at 
the "optimal phonon" emission energyi^ This energy is 
on the order of hs/Xd, which is the energy of an acoustic 
phonon with the wavelength equal to the dot size A^- 

The Coulomb interaction between the electrons has a 
number of effects on the spin relaxation. Firstly, the 
Coulomb interaction affects the relaxation rates via the 
singlet-triplet energy Ets, because the relaxation rates 
depend strongly on Ets- The singlet-triplet energy Ets 
is sensitive to the Coulomb interaction and it is gener- 
ally smaller than in the non-interacting case. Moreover, 
the singlet-triplet crossing {Ets = 0) is possible in finite 
orbital magnetic fields only if the Coulomb interaction 
is present. Secondly, the Coulomb interaction modifies 
the form-factor of the spin-phonon interaction and in- 
troduces a new frequency scale (for strong Coulomb in- 
teraction) in the relaxation rates, uj s/{r), where (r) 
is the average distance between the electrons and s is 
the speed of sound. Finally, the spin admixture mecha- 
nism that allows for the singlet-triplet relaxation to occur 
is essentially different in the limits of weak and strong 
Coulomb interaction. For weak Coulomb interaction, su- 
perimposing states inside the singlet-triplet subspace is 
sufficient for obtaining efficient singlet-triplet relaxation. 
For strong Coulomb interaction, this mechanism becomes 
inefficient. Despite the fact the states in the singlet- 
triplet subspace are strongly admixed to each other, the 
orbital parts of the wave functions of the singlet and 
triplet become similar to each other. More precisely, the 
electron charge distributions in the singlet and triplet 
states converge onto each other in the limit of strong 
Coulomb interaction. Therefore, a different admixture 
mechanism, which was neglected in this paper, is respon- 
sible for singlet-triplet relaxation in the limit of strong 
Coulomb interaction. 

We modeled the quantum dot by a circularly symmet- 
ric harmonic confining potential and took into account 
the Coulomb interaction with a high accuracy using a 
variational method. The quantum dots used in experi- 
ment might, of course, not have a harmonic confining po- 
tential and not even have circular or inversion symmetry. 
We, therefore, indicate here which of our results change 
if the quantum dot confinement assumes a different form. 
Our results about the spin relaxation rates remain quali- 
tatively the same for any confining potential. This refers 
to the non-monotonic dependence of the spin relaxation 
rates on the magnetic field shown in Figs. [TO] and [TT] 
The expressions in Eqs. PU^ . piI5)l and piH)) . as weU 



as in Eq. (|138p . are valid only for circularly symmetric 
confining potentials. If the confining potential has no cir- 
cular symmetry, then an additional integration over the 
in-plane angle </> of the emitted phonon is present in each 
of these expressions. The dependence on the form of the 
dot confinement enters in the relaxation rates through 
the form-factors TsriQw) and TT{q\\)- In Eqs. (|98p and 
(|136l) . the form- factors were evaluated for the harmonic 
confinement and in the presence of the Coulomb interac- 
tion. An important difference between confining poten- 
tials with and without center of inversion for the singlet- 
triplet relaxation was mentioned below Eq. (|107p . In 
the limit q\\ 0, one obtains ^sriQw) ^ 9^ for con- 
fining potentials with a center of inversion. If no cen- 
ter of inversion is present, then ^ST^Qx'^Qy') ~ Qx',y'^ 
i.e. the singlet-triplet form-factor is proportional to the 
first power of As a result, in the limit of vanish- 
ing phonon transition energy, w — > 0, the singlet triplet 
relaxation rates for the case of a confinement without 
inversion symmetry have a weaker w-dependence (two 
powers of w less) than for the case of a confinement with 
inversion symmetry. The scale of w where the crossover 
takes place is determined by the degree of asymmetry of 
the potential. The dependence of the avoided splitting 
gaps A± on the direction of the magnetic field is also 
sensitive to the form of the confining potential. Two lim- 
iting cases, namely of a circularly symmetric dot and of 
an elongated dot (or a double dot) , have been considered 
in Sec. IVBl (cf. Fig. H and Fig. The matrix elements 
of the spin-orbit interaction have general forms given in 
Eqs. dill), ((251), (HH), and §0^, which are vahd for an ar- 
bitrary shape of the dot confinement. They can be used 
to express A± in terms of the matrix elements (n|^''|n') 
of a dot with an arbitrary confining potential, as long as 
the characteristic extension of the orbital wave function 
is much smaller than Xso- In Sec. IV A[ we have not as- 
sumed any specific form of the quantum dot potential, 
however, we restricted our analysis to quantum dots that 
have no orbital degeneracies of the two-particle levels. 

Another important assumption made in this paper is 
that the spin orbit interaction is linear in the electron 
momentum. This type of the spin-orbit interaction is 
most commonly used for GaAs quantum dots, since the 
confinement in the normal to the 2DEG direction {z- 
axis) is strongest. In Sec. IV Al we showed that, for this 
type of spin-orbit interaction at the leading order, no 
spin relaxation occurs within the singlet-triplet subspace 
without the presence of a Zeeman splitting. As a re- 
sult, the lowest-order singlet-triplet relaxation rates are 
proportional to the second power of the Zeeman energy, 
in agreement with the experiments on triplet-to-singlet 
relaxation in quantum dotsi^ 

In conclusion, we found strong variations of the spin- 
relaxation rates (by many orders of magnitude) across 
the singlet-triplet crossing, with extremely large values 
(1 — 100 s) of spin lifetimes at the avoided crossing points. 
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APPENDIX A: VARIATIONAL PARAMETERS u 
AND 7 

Here, we study the variational parameters uj and 7 
that minimize the energy egm in Eq- (|13[) . We derive 
asymptotic expressions for uj and 7 in the limits of strong 
{\/a*g ^ 1) and weak (A/a|j ^ 1) Coulomb interaction. 
Equating to zero the derivatives of epm in Eq. (jl3p with 
respect to lu and 7, we obtain two equations 



m 

V 



X r(l/2 + j;) 



2a% 



xr{x) 



0, (Al) 



2 4 
X y - 



+ [vl/(l + - $(1/2 + x)] = 0, 

(A2) 

where i = a/ + 7, y = X/X ~ y/uj/oj, and ^'(a;) is 
the digamma function. Solving Eqs. (lAip and (|A2p with 
respect to x and ?; gives the variational parameters 7 and 
u as functions of the Coulomb interaction strength X/ a*g. 

We consider first the case m = 0. For weak Coulomb 
interaction (A/a|j <C 1), we expand Eq. (jA2p in terms of 
a; <§; 1 and obtain 



2 In 2 



41n 2 + 7r2/3 



(A3) 



Note that, for X/a*^ < j/^/20rln2, Eq. ([XS]) gives a; < 0, 
whereas by definition x = ^ > 0. In this case, the 
minimum of energy in Eq. (|13p is achieved at 7 = and 
thus, Eq. (jA2p should be replaced by a; = for X/a^ < 
y^ /2^\n2. Setting a; = in Eq. (|Aip . we obtain an 
equation for y, 



-2/- 1 = 0, 



(A4) 



which is vahd for X/a*g < y^ /2y^\n2. Consider- 
ing simultaneously the equality X/a*g — ?/'^(2^/7r In 2)~^ 
and Eq. (|A4p . wc find that X/a^ has a critical value, 
(^/«s)c = L, given by 
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0.57, 



(A5) 



at which a; and 7 are non- analytic functions of X/a*g. 
Clearly, such a critical point is not present in the ex- 
act eigenstates of the Hamiltonian ([9]), and is an arti- 
fact of the variational ansatz we use. In the interval 
< X/a*g < {X/a*g)c, we have 7 = and lu = Lu/y'^, 
where y is the positive solution of Eq. (|A4p . Although 



we can solve Eq. (jA4p analytically for y, it is more con- 
venient to present here an expansion for uj in terms of 
X/a%, 
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(A6) 



Note that, since X/a*^ < {X/a*g)c ~ 0.57, Eq. (jA6p con- 
verges with a good accuracy; e.g., at X/a*^ — {X/a*g)c the 
four terms in the series (|A6p give uj/lu Ri 0.797, whereas 
the exact solution reads Cj/uj ~ y/l — 1/4 In 2 « 0.799. 

To the right of the critical point (jASj) . both Eqs. (jAT|) 
and (IA2D are valid for m = 0, and we can use them to 



study two limiting cases: (i) the neighborhood of the 
critical point (jA5p for X/a^ > {X/a*g)c and (ii) the limit 
of strong Coulomb interaction {X/a*g ^ 1). Excluding 
X/a*g from Eqs. (|Aip and (|A2p . we express ?/ as a function 
of X only. 



2[*(l + a;)-«'(l/2-(-a;)]' 



(A7) 



Next we substitute y from Eq. (|A7p into Eq. (|A2p and 
obtain 

f = ^ii(i±^[*(l + .)-*(l/2 + ..)f^^ 
A 1 (1 + xj 

X {(1 + x) [*(1 + a;) - *(l/2 + a;)] - 1/2}'^^ . (A8) 



Equation jASj) gives a*g/X as a function of a;. The in- 
verse function allows one to find a; as a function of 
a|j/A, and then, from Eq. (|A7p . also y as a function 
of a^/A. Expanding Eqs. and fXT)) for a: < 1, 

which corresponds to the case (i), we find for ^y/7 = x 
and oj/oj = l/y^ in leading order, 
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(A9) 
(AlO) 



where is given in Eq. (|A5p . {uj/lj)c = \/l — 1/4 In 2, 
and the coefficients Co and Ci read 



(2/^)i/"(ln2)3/4(41n2- 1)1/4 
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0.014. 



(A12) 



Note that the value Cj/uj = {Cj/uj)c at A/a|j = is the 
minimum of Cj/lo &s a. function of A/a|j [see also Eq. (|A6P 
and the text below it] . 

Next we turn to the case of strong Coulomb interaction 
{X/a*g » 1). Here, we have a; 1 and by expanding 
Eq. (jA7p in terms of 1/a; we obtain 
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8j) for large x, and inverting the ob- 

2 



Expanding Eq. 

tained series by means of iteration, wc find for j = x 
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Then, Eq. (jA13|) can be rewritten as follows 
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To summarize for the case m — 0, the variational pa- 
rameters 7 and a) as functions of the Coulomb interac- 
tion strength A/a|j are given, respectively, by 7 = and 
Eq. in the interval < A/a^ < ^c, by Eqs. and 
(|A10|) in the neighborhood of A/a|j = (for X/a% > ^c), 
and by Eqs. (|A14|) and (jAlSp in the hmit X/a% > 1. 

We consider now the case \m\ > 1, and without loss 
of generality we assume m — \m\. Here, we find that uj 
and 7 are analytic functions of \/a*g- Wc thus consider 
only the limit of weak and strong Coulomb interaction. 
Proceeding similarly to the previous case, we find from 



Eqs. 
and 



T]) and (|X2|l the following relations [cf. Eqs. 
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where A5'(x) = ^'(l-fx) — '^{1/2 + x) and x takes values 
in the interval to < a; < 00. 

In the limit of weak Coulomb interaction (A/a|j <C 1), 
we expand the right-hand side of Eqs. (|A16|) and (|A17|) in 
terms of x — to « ^ /2m <IC 1, and obtain in leading order 
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In particular for 

5/4)(A/a^3) 



where A^{m) = 2 In 2 - EfeLi 
the case |to| = 1, we obtain 7 = •y7r(21n2 
and Lj/uj = 1 - 0F(3/2 - 21n2)(A/4a|;). 

In the limit of of strong Coulomb interaction (A/a|j ^ 
1), we expand the right-hand side of Eqs. (|A16|) and (|A17|) 
in terms of a; 3> 1, and as a result we obtain 
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FIG. 12: (Color online) (a) The variational parameters 7 as 
a functions of the Coulomb interaction strength X/a'g for the 
singlet m = Q and the triplet \m\ = 1 states, (b) The ratio 
iD/oj as a function of the decimal logarithm of A/a^ for m — 
and |m| = 1. To find 7 and uj/uj we minimize the ground state 
energy in Ea.(|13p. see Appendix 1X1 for details. 



interaction, the radial part of the wave function in Eq. ([T]) 
[see also Eq. (jl2p ] is weakly depending on the quantum 
number to. 

In Fig. [13 we plot 7 and uj/uj as functions of A/a|j 
and logio(A/a^), respectively, both for the singlet state 
m — and the triplet states to = ±1. Note that, for 
the singlet state to = 0, the parameters 7 and uj behave 
non-analytically at the point A/a|j = ~ 0.57, where 
is given in Eq. (|A5p . Nonetheless, the energy of the 
TO = state, calculated according to Eq. (fT5|) . is a smooth 
function of A/a|j, i.e. it has a continuous first derivative 
with respect to A/aJj. 



APPENDIX B: 



CROSSOVER FORMULAE FOR 5 

AND ujZ 



The singlet-triplet splitting Ets defined in Eq. 
contains an interaction parameter 

^ _ £0,-1 ~ gpo 

huj ' 



(Bl) 
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where the energies e„m are the eigenvalues of T-Lm in 
Eq. ([5]). The parameter 5 shows by how much the en- 
ergy levels are renormalized by the Coulomb interaction. 
In Fig. [infa), we plot the parameter (5 as a function of 
A/a|j, calculated with the help of Eq. (fT3|) (solid curve). 
For analytic calculations, it is convenient to use the fol- 
(A2b)wing formula 
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(B2) 



Note that in Eqs. (|A20[) and (|A2ip the dependence on to 
arises only in the last terms. Thus, for strong Coulomb 



where h 
We obtained Eq. 
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FIG. 13; (Color online) (a) The parameter 5 as a function 
of the Coulomb interaction strength X/a*g. The dotted curve 
shows S calculated according to Eq. (|B2|I . (b) The singlet- 
triplet degeneracy occurs at a megnetic field B* {Ets = 0). 
This figure shows the corresponding lo'/loq as a function of 

and Ao = ^2h/m*LUQ. The 

(El. 



Xo/ag, where cu* — eBl/m'c 

dotted curve corresponds to the fitting function in Eq. 



for the cases of strong and weak Coulomb interaction. 
Therefore, Eq. (|B2[) is exact in the hmits of weak and 
strong Coulomb interaction, and is accurate within 7% 
in the crossover region {h^^ < \/a*g < c^^/"^). We plot 
5{\/a*B) given by Eq. in ¥ig.^a) (dotted curve). 

The singlet-triplet splitting Ets in Eq. ([H]) goes to 
zero at a value of cyclotron frequency uj*^. In Fig. [TST b). 
we plot the ratio ijJ*/u!o as a function of a|j/Ao; the dot- 
ted curve corresponds to using the following crossover 
function 



ri/3 



- 1 HX^/al) 



where T{x) = (l -h O.llx^) / (l + O.Sx^/^) . The 
crossover function (|B3p was obtained, similarly to 
Eq. (|B2p . with the help of asymptotic expressions for Cj 
and 7 from Appendix \K\ for the limits of weak and strong 
Coulomb interaction. 



APPENDIX C: FIDELITY OF VARIATIONAL 
METHOD 

We discuss now the accuracy of our variational 
method. We rewrite Eq. ^ in the following form 



where Hm has the eigenvalues £„, 



(CI) 



\/ rr? -I- 7) and the eigenfunctions given in Eq. (fT2|) . The 
perturbation 
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FIG. 14: (Color online) Infidelity 1 — F of the ground state 
of Tim for |m| = 0, 1. We used Eq. (IC5[l with the matrix 
elements of V evaluated in appendix [D] Note that both axes 
are logarithmic (with decimal base). 



is expected to be "small" when the variational method 
works, with "small" meaning that the off-diagonal ma- 
trix elements of V are much smaller than the level spacing 
(the diagonal part of V needs not be small). We can show 
that indeed V <C Tim in the case A/a|j <C 1, and also in 
the case A/a|j 3> 1 for transitions involving the ground 
state. As for intermediate strengths of the Coulomb in- 
teraction {X/a*g ^ 1), there is no small parameter in the 
problem, which can justify the assumption V <^ Ti-m- 
Nonetheless, we calculate the matrix elements of V (see 
Appendix[D| between the ground and excited states, and 
find that they are numerically small, as compared to the 
level spacing 2huj. In particular, we find that the matrix 
element {fim\V\fom) is identically zero, if the variational 
parameter lu minimizes the energy of the ground state 
in Eq. (|13p [see also Eq. (|Aip ]. Furthermore, we deter- 
mine an upper bound on the absolute value of any of 
the remaining matrix elements (n > 2). We find that, 
in both m = and |m| = 1 cases, the matrix element 
{f2.m\V\fom) has the largest value in the crossover re- 
gion X/a'g ~ 1. If measured in terms of the level spacing 
2huj, this upper bound is approximately 0.1 and 0.03, for 
the m = and |m| = 1 cases, respectively. Next, we 
characterize the accuracy of our variational method in 
terms of the fidelity 



F 



(C3) 



where 'ii{r) = fomir) is the exact ground state of Hm, 
and 'l'(r) = /om(^) is the corresponding wave func- 
tion calculated with the help of the variational method. 
We estimate the fidelity, using the perturbation theory 
expansior 



fomir) = fomir) + %^^5%^/„™(r) + ... , (C4) 
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where we retain terms up to the second order in V (not 
shown). Thus, we obtain the infidehty 



Finally, for the term V" = jm*{uj'^ — (I'^)r^, we obtain 



1 - F = 



^ OO 

n=2 



{fnm\V\f 071 



2nhib 



(C5) 
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* \2 



(D4) 



where we dropped the n = 1 term, due to its zero con- 
tribution. The matrix elements Vnn' = {fmn\V\.fn'm) are 
calculated in Appendix [Dl In Fig.[T31 we plot lg(l — F) as 
a function of lg(A/a|5) for \m\ — 0, 1, with lg(a;) being the 
decimal logarithm of x. Figure [14] shows that our vari- 
ational method is fairly accurate, even in the crossover 
region X/a*g ~ 1. Finally, we note that in Eqs. (IC4p and 
(|C5|) we did not include the diagonal part of V in the en- 
ergy denominators. More accurately, one has to take this 
part into account as well, since Vnn can be comparable 
to the level spacing huj. However, it turns out that the 
difference Vnn — Voo > monotonically increases with n 
in the whole range of X/ag. Therefore, when we replaced 
2nhw 2nhu) + Vnn — Vqq in Eq. (jCSp . we obtained only 
an insignificant reduction of 1 
scale of Fig. [T3] (not shown) . 



F for X/a*g ~ 1 on the 



APPENDIX D: MATRIX ELEMENTS OF V 

We calculate first the matrix elements {fnm\V\fom) 
(n > 0), with V given in Eq. (|C2p and the wave func- 
tions fnm {r) in Eq. (fT2|) . We divide the perturbation into 
three terms, 

Ko = (/nrnl^l/o™) = <0 " + ^nd, (Dl) 

corresponding, respectively, to the terms proportional 
to l/r, l/r^ and in Eq. ((C2)) . For the term V = 
fi? /m*a*gr, we obtain 



where 5nn' is the Kronecker J-symbol. We note that 
in the case rt = 1, we have (/im|^|/om) = 0, due to 
Eq. (IHIl . 



Next, we consider the general case, 



Vnn' = (/„™ |y = " K:'n' + Kn' , (D5) 

with the same division of the perturbation into three as 
in Eq. (jPip . For the first part, we obtain 
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where (x)„ — r(a; -I- n)/r(a;) is the Pochhammcr symbol. 
For the second part, we obtain 
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(D7) 

We note that for n' = n we have Vnn — ^'^^ all values 
of n. Finally, for the last part, we obtain 



ft2 (2n- l)!!2-"r(i- 1/2) 
m*a*sX y^n\T{t)T{n + 1) 



(D2) 



Here and below we use the notation t — 1 + \J + 7. 
For the term V" — fi,^7/m*r^, we obtain 



V , = 



m*X 



* \2 



{up — Up') {2n + t)Snn' ~ \/ n'[n 4- 1) Sn 



'y/n{n' + t) 5n,n' + l 



(D8) 



V" - 



*A2 (i-1) V r(n + i)' 



We note that for n' ^ Eqs. (lD6|, jOTl and (|D8]) coin- 
(D3) cide with Eqs. JDS]), (|D3]) and (lD4|, respectively 
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